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We present a high statistics study of the light hadron spectrum and quark masses in QCD with 
two flavors of dynamical quarks. Numerical simulations are carried out using the plaquette gauge 
&J}, action and the 0(a)-improved Wilson quark action at /3 = 5.2, where the lattice spacing is found to 

^ < be a = 0.0887(11) fm from p meson mass, on a 20^ x 48 lattice. At each of five sea quark masses 

' corresponding to mps/my — 0.8-0.6, we generate 12000 trajectories using the symmetrically pre- 

conditioned Hybrid Monte Carlo algorithm. Finite spatial volume effects are investigated employing 
12"^ X 48, 16'' X 48 lattices. We also perform a set of simulations in quenched QCD with the same 
lattice actions at a similar lattice spacing to those for the full QCD runs. In the meson sector we 
find clear evidence of sea quark effects. The J parameter increases for lighter sea quark masses, and 
^K^J ■ the full QCD meson masses are systematically closer to experiment than in quenched QCD. Careful 

JiT ' finite-size studies are made to ascertain that these are not due to finite-size effects. Evidence of sea 

' quark effects is less clear in the baryon sector due to larger finite-size effects. We also calculate light 

Q ! quark masses and find m^(2 GeV) = 3.223 (to.oM) ^eV and m^{2 GeV) = 84.5 MeV 

' which are about 20 % smaller than in quenched QCD. 
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I. INTRODUCTION 



^ il Lattice QCD calculation of the light hadron mass spectrum has witnessed significant progress in recent years 

^ ■ H S H. In the quenched approximation in which the quark vacuum polarization effects are ignored, the 
' CP-PACS Collaboration performed a precise calculation of hadron masses, in which the estimated accuracy 
reached the level of a few percent in the continuum limit 4]. They found that the quenched spectrum shows 
a significant and systematic deviation from experiment; the K*-K hyperfine splitting is smaller by about 10% 
^ . , than experiment. The decuplet baryon mass splittings are also small, and the octet baryon masses are themselves 
' smaller than experiment. 

Since this work, the focus of efforts has shifted toward full QCD simulations including vacuum polariza- 
tion effects of dynamical quarks. A number of simulations now exist, pursued by the SESAM-T^L |1 
UKQCD CP-PACS |l3,|i3 and QCDSF-UKQCD Q Collaborations for two flavors using the Wilson-type 

quark action, and by the MILC Collaboration [TsLll^ for two and three flavors using the Kogut-Susskind(KS) 
quark action. In particular the CP-PACS calculation [loL HH made a first attempt toward execution of the 
chiral and continuum extrapolations within the single set of simulations, as pioneered by the GFll Collabora- 
tion ^5 in their quenched spectrum study. The chief finding of this work was that the K*-K hyperfine splitting 
agrees much better with experiment in two-fiavor full QCD than in quenched QCD, and that light quark masses 
decrease by about 25% by the inclusion of dynamical u and d quarks. 

A subtle point with the CP-PACS results is that the dynamical sea quark effects become manifest only after 
the continuum extrapolation. Further studies are required to consolidate effects o f dy namical sea quarks. The 
CP-PACS simulation used a renormalization group (RG) improved gauge action llq, but the 0(a)-improved 
Wilson quark action 17] with only tadpole-improved jlgj clover coefficient. This leaves scaling violation of 
0{g^a). The use of non-perturbatively determined clover coefficient removing all of 0{a) errors should be much 
better to control the continuum extrapolation [T9|. Studies along this direction were previously carried out by 
the UKQCD and QCDSF Collaborations using the plaquette gauge action. However, sea quark effects are not 
clear in their results of hadron masses, albeit encouraging evidence is seen in the static quark potential |9l.ll2||. 

In the present work, we explore sea quark effects in hadron and quark masses in two-flavor QCD using the 
plaquette gauge action and the non-perturbatively 0(a)-improved Wilson quark action. Our simulations are 
performed at a single lattice spacing a^^ ~ 2 GeV at /3 = 5.2 using a 20^^ x 48 lattice. We also carry out 
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calculations in quenched QCD with the same action and similar simulation parameters to those in full QCD in 
order to make a direct comparison between full and quenched QCD. Preliminary results of these calculations 
have been reported in Refs. [20ll2ll23. |2^ . 

We pay particular attention to two points which are important for an unambiguous identification of sea quark 
effects. One is the finite spatial volume effect whose magnitude is believed to be more pronounced in full QCD 
simulations than in quenched QCD (24i | . An increase of hadron masses due to this effect could mimic sea quark 
effects. There are only a few studies of finite-size effects in full QCD for the Wilson- type quark action [aH]- 
This leads us to perform a systematic investigation of finite-size effects employing 12"^ x 48, 16'^ x 48 and 20^ x 48 
lattices. 

Another point is the chiral extrapolation. With currently available computer power and simulation algorithms, 
the sea quark mass which can be explored with the Wilson-type quark action is limited to values corresponding 
to 77ips,sca/'7iv,soa ^ 0.6. The loug extrapolation to the physical u and d quark masses may involve sizable 
systematic errors, potentially blurring, or artificially enhancing, sea quark effects. This can be avoided if one 
examines sea quark effects at the quark masses actually simulated. In this study, therefore, sea quark effects 
are examined in detail not only at the physical quark mass but also at our simulation points. 

We have also made efforts to accumulate high statistics of 12000 trajectories each at five values of sea quark 
masses. Our implementation of the symmetric preconditioning of the lattice clover-Dirac operator | i25i . i2^ 
speeded up the configuration generation by a factor two by allowing a doubly larger step size over the even-odd 
preconditioning . 

This paper is organized as follows. We describe details of configuration generation in full and quenched QCD 
in Sec. [H] Method of measurement of hadron masses and the static quark potential is explained in Sec. Illll 
Finite-size effects on hadron masses are discussed in Sec. lIVI Section^is devoted to detailed description of the 
chiral extrapolation of our hadron mass data. We examine sea quark effects in light hadron masses in Sec. IVII 
Results of the decay constants and quark masses are presented in Sees. IVIll and IVIlTl Our conclusion is given 
in Sec.HXI 



II. SIMULATION METHOD 



A. Simulation parameters and algorithm 



We carry out numerical simulations of lattice QCD with two flavors of degenerate dynamical quarks which 
are identified with the up and down quarks. We use the standard plaquette action for gauge fields defined by 

-5^ = ^ E (1) 
where Ux,fj,i^ is the product of gauge link variables Ux.fi around the plaquette given by 

Ux,iiv — Ux,fj.Ux+fi,iyUl^^^^Ul^^. (2) 

The 0(a)-improved Wilson action defined by 

Sq = ^qxDxyqy (3) 

x,y 

Dxy = S^y - Ky^ {{l - 7^) Ux,t^Sx+f^,y + (1 + 7^) Ux.n^x^y+p.} - -^Kcswcr^uFx^^^Sxy (4) 

is used for the quark part. The field strength tensor on the lattice is defined by 

1 

8i 



where (h.c.) denotes the hermitian conjugate of the preceding bracket, and ct^^ ~ (*/2) [7u, 7J. The clover 
coefficient csvv is set to the non-perturbative value determined by the ALPHA Collaboration [27l |. 

Our simulations are performed at a single value of (3 — ^.2. The lattice spacing fixed from nip at the physical 
sea quark mass is found to be 0.0887(11) fm. Our value of (3 is slightly off the range /3 = 12.0-5.4 where the 
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ALPHA Collaboration carried out a non-perturbative determination of csw- We set csw = 2.02 by extrapolating 
their parametrization formula of the non-perturbative csw as a function of the bare coupling. We performed 
an independent non-perturbative determination of csw at (3 = 5.2 and confirmed that our preliminary result 
Csw = 1-98(7) is consistent with 2.02 within the error p8|. 

We employ three lattice sizes that differ in spatial volumes, A^^ x A^t = 12^ x 48, 16^ x 48 and 20^ x 48. The 
hadron spectrum and quark masses at the physical point are calculated using the data on the largest lattice. 
The data on the two smaller lattices are used to investigate finite-size effects on hadron masses. 

On each lattice size, we adopt five values of the sea quark mass corresponding to the hopping parameter 
J^sca = 0.1340, 0.1343, 0.1346, 0.1350 and 0.1355. This choice covers tops, sea/ wy, sea = 0.6-0.8, and enables 
us to extrapolate our data to the physical sea quark mass. These simulation parameters are summarized in 
Table D 

We note that the UKQCD Collaboration also performed a set of simulations using the same lattice action at 
a~^o±2 GeV 0. There are, however, some differences in the choice of /? and A'sca^ The UKQCD simulations 
shift /9 with the sea quark mass keeping the Sommer scale ro/a |29j | fixed, while our simulations are performed at 
fixed f3. Another difference is the range of the sea quark mass covered in the two simulations. We explore light 
sea quark masses down to mps.soa/TOv.soa — 0.6, whilst UKQCD's lightest point is around mps.soa/TOv.sca — 0.7. 
Although the UKQCD Collaboration made another simulation at a smaller sea quark mass tops, soa/TOy.soa — 0.6 
at a spatial extent of A^^ ~ 16 (A^sO— 1-6 fm), finite-size effects seem to be significant there (see discussion in 
SecHg. 

Gauge configurations are generated using the Hybrid Monte Carlo (HMC) algorithm [s^,!^- We use simu- 
lation programs with three variants of HMC for the 0(a)-improved Wilson action: 

• HMC with the even/odd preconditioning |3^ only for the inversion of the quark matrix D^y This 
algorithm is used in the simulations on the 16'^ x 48 lattice. 

• HMC with the asymmetric preconditioning for the lattice action (A-HMC) [25ll2^l33| . Whole simulations 
on the 12"^ x 48 lattice are performed with this algorithm. 

• HMC with the symmetric preconditioning for the action (S-HMC) |2^ , which shows the best perfor- 
mance among the three algorithms. 

Our main simulation on the 20"^ x 48 lattices is initially started with the A-HMC algorithm, but is later switched 
to the S-HMC to speed up the calculations. The trajectory length in each HMC step is fixed to the unit length. 
We use the conventional leap-frog integration scheme for the molecular dynamics equation. The step size At is 
chosen to achieve an acceptance of 60-80%. 

The even/odd preconditioned BiCGStab algorithm [s^l is used for the quark matrix inversion to solve the 
equation D^yGy — B^- We take the stopping condition of the form \\DG — B\\ < A in the HMC program. 
A modified form \\DG - B||/||B|| < A is used in the A-HMC and S-HMC programs. The value of A in the 
evaluation of the fermionic force is determined so that the reversibility over unit length is satisfied to a relative 
level better than 10"^'^ for the Hamiltonian. We use a stricter stopping condition in the calculation of the 
Hamiltonian in the Metropolis accept /reject test. Table ^ shows our choice of A together with the average 
number of the BiCGStab iteration in the quark matrix inversion for the force calculation, A^inv 

We accumulate 12000 HMC trajectories at each sea quark mass on the 20"^ x 48 lattice. The statistics on 
smaller lattices are 3000 trajectories. Measurements of light hadron masses and the static quark potential are 
carried out at every 10 HMC trajectories. Details of the measurement method will be described in the next 
section. 

All simulations are performed on the parallel computer HITACHI SR8000 model Fl installed at KEK. This 
machine consists of 100 nodes and has a peak speed of 1.2 TFLOPS and 448 GB of main memory in total. 
The CPU time needed per unit HMC trajectory on the full machine is listed in Table ^ The total time for 
configuration generation on each lattice size is 8.6 days on 12"^ x 48, 58 days on 16"^ x 48, and 130 days on 20'^ x 48 
lattices. Additional 100 days are spent for the measurement of the hadron masses and the static potential. 
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B. Simulation in quenched QCD 



While many calculations of the hadron spectrum have been performed in quenched QCD, comparisons between 
our full QCD results and quenched results from other simulations may be subject to systematic uncertainties 
due to the difference in the simulation details. We therefore carry out a set of quenched calculations of the 
hadron spectrum using the same lattice actions and simulation parameters as those for full QCD runs. 

Our simulations are performed at /3 = 6.0, where the lattice spacing fixed from nip equals 0.1074(14) fm. We 
take csw = 1.769 which is the value determined non-perturbatively by the ALPHA Collaboration js^- Three 
lattice sizes 12^ x 48, 16"^ x 48 and 20"^ x 48 are employed in order to investigate finite-size effects. 

Gauge configurations are generated with a combination of the heat-bath and the over-relaxation algorithms. 
We call four heat-bath sweeps with a succeeding over-relaxation step an iteration. We accumulate statistics 
of 60000 iterations on each lattice size. Hadron masses and the static potential are calculated at every 200 
iterations. 



III. MEASUREMENT 



A. Hadron Masses 



In measurements in full QCD, we use six values of the valence quark mass corresponding to the hopping 
parameter K^^i^i{i = 1,...6) = 0.1340, 0.1343, 0.1346, 0.1350, 0.1355 and 0.1358, which cover the range of 
'TipS,vai/TOv,vai — 0.5-0.8. At cach sea quark mass, therefore, there is one value of i^vai,j, which equals Kgca 
and is identified as the light quark mass. Other five values of K^ai^i correspond to the mass of strange quarks 
treated in the quenched approximation. In the following, we use the abbreviation "diagonal data" to represent 
hadron correlators or masses with a quark mass combination in which all valence quark masses are equal to the 
sea quark mass. 

We employ meson operators defined by 

M{x) = qi^Wf , r = /, 75, 7m, 757m: (6) 

where / and g are flavor indices and x is the coordinates on the lattice. Meson correlators {M{x)M{Oy) are 
calculated for the following eleven combinations of valence quark masses 

(i4rval,i,i^val,i) (i = 1, ■ • -6), 

(7) 

{Ksca, K^&u) (l = 1, . . . 6, Ksca K^&u)- 

The former is identified with a degenerate light or strange meson and the latter with a non-degenerate light- 
strange meson. This choice of the valence quark masses enables us to calculate the full spectrum of strange and 
non-strange mesons. 

For baryons, we use the same operators as those employed in Ref. [ll| . Namely, the octet baryon operator is 
defined as 

0/9'^(x) = e-"- (qy^''^C-i,qi^'^') qi'^'>'^, (8) 

where a, 6, c are color indices and C = 7472 is the charge conjugation matrix. We measure baryon correlators 
with two types of flavor structure (S- and A-like baryons), 

S : --L^ol/^ls + ols"]/), (9) 

A : -L (^0[/'^l9 - 0^3h]f _ 2o[/9]''j , (10) 

where Ol/sl'* = O^s'' — O^^^. Decuplet baryon correlators are calculated using an operator defined by 

= ,abc (g(/)-^C7M9i^"') gi"'^ (11) 
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with symmetrized flavor structure 

^ (^Jjfgh £)hfg _|_ jjghf _|_ j^fhg _|_ jjgfh ^ £)hgf\^ 



(12) 



We take quark mass combinations of {K^a\,i,K^s\A,K^aia) (« = l,---6), (i^sca, -fi^vai,j, -fi^vai.i) and 
(/Csca, i^sca, i^vai/i) (i = 1, . . . 6, i^soa 7^ -f^vai,i) for the baryou correlators. 

In order to construct the smeared hadron operators, we measure the wave function of the pseudoscalar (PS) 
meson 



Ex(g(x,075g(x + r,t)P(O,0)t) 
E.(P(x,i)P(O,0)t) 



(13) 



where P is the PS meson operator, Eq. ©, with P = 75 and t fixed to 12. The measurement is performed 
at each sea quark mass and lattice size using a subset of gauge configurations (30 configurations every 100 
trajectories). We parameterize (j){v) using a polynomial approximation 0(r) = 1 + X]n=i s'^"!''!" ^'^^ 
the smearing function. We employ three types of the meson operator : i) local operator, ii) smeared operator 
M{x) = 0(r)g(x, i)Pg(x + r, i), in) doubly smeared operator M{x) = '^^^ ^, (p{r)(p{r')q{x + r, t)rq{x + r', t). 
Additionally, we use "triply smeared operator" 

O^s^x) = 0(ri)0(r2)</>(r3)e"''^(g"^(x + ri,i)CT<?''(x + r2,i))9^(x + r3,i) (14) 

ri,r2,r3 

for baryons. Hadron correlators are measured with a) point source and sink operators, b) smeared source and 
point sink, and c) smeared source and sink operators. We fix configurations to the Coulomb gauge, since b) 
and c) are not gauge invariant. 

We observe that, when valence quarks are lighter than sea quarks, the hadron correlator takes an excep- 
tionally large value on a small number of configurations. This might be caused by a fluctuation of the lowest 
eigenvalue of the Dirac operator of the 0(a)-improved Wilson action. If the PS meson correlator on the z-th 
gauge configuration takes a value larger than 20 times the statistical average, which is evaluated without that 
configuration, at a certain time slice 

N r 
on "coat 

(P(a:)P(0)t), > iPi^)m^)k (15) 

where iVconf is the total number of configurations, we consider it as an exceptional configuration and remove it 
from the following analysis. The number of the removed configurations is given in Tabic ^ 

In order to reduce the statistical fiuctuation of hadron correlators on the 20"^ x 48 lattice, we repeat the 
measurement for two choices of the location of the hadron source, t^i-c — 1 and A'^t/2 + l(= 25), and take the 
average over the two sources: 

i ((M(i,,.e + t)Af(t,,e)^),^^^^i + {M{t,rc + t) M (U,,)^ ) ^^^^^^'^ . (16) 

We find that this procedure reduces the statistical error of hadron correlators by typically 20%, which suggests 
that the statistics is increased effectively by a factor of 1.5. For further reduction of the statistical fluctuation, 
we take the average over three polarization states for vector mesons, two spin states for octet baryons and four 
spin states for decuplet baryons. 

Figures ^-01 show examples of effective mass plots. We find that the best plateau of the effective mass is 
obtained from hadron correlators with the point sink and the doubly smeared source for mesons and the triply 
smeared one for baryons. Therefore, hadron masses are extracted from these types of correlators. 

We carry out fits to hadron correlators by taking account of correlations among different time slices. A 
single hyperbolic cosine form is assumed for mesons, and a single exponential form for baryons. The lower cut 
of the fit range tmin is determined by inspecting stability of the fitted mass. The upper cut (tmax) dependence of 
the fit results is small and, therefore, we fix i^ax to Nt/2 for all hadrons. Our choice of fit ranges and resulting 
hadron masses are summarized in Tables IXXIXhIXXXIVI in Appendix A. Statistical errors of hadron masses 
are estimated with the jack-knife procedure. We adopt the bin size of 100 trajectories by inspecting the bin size 
dependence of the jack-knife error as discussed in Sec lIII Cl 
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In Fig.|31 we test double exponential fits to extract hadron masses at iiTgca = ^vai = 0-1355. While these fits 
are unstable and lead to a large error for the mass of the first excited state, the result for the ground state 
mass is consistent with that from the single exponential fit. The situation is similar at other sea and valence 
quark masses. This suggests that the hadron masses in Tables IXXIXI - IXXXIVI and the light hadron spectrum 
calculated from these results have small contamination from excited states. 

Hadron correlators in quenched QCD are calculated in an analogous manner. We use six values of ifvai, 
0.13260, 0.13290, 0.13331, 0.13384, 0.13432 and 0.13465, corresponding to mps,vai/mv,vai=i 0.50 -0.80 and the 
hadron operators and smearing procedure same as those in the full QCD study. A difference is that we can take 
more combinations of valence quark masses than in full QCD, since any value of the six valence quark masses 
can be identified with either light or strange quark mass. We take all combinations (i^vai.i, J^^vai.j) {hj — 1, ■•, 6) 
for mesons, and somewhat restricted choices (-ftTvai,!, ^vai,i, ^vaij) {hj = l,--,6) for baryons. Statistical errors 
are estimated with the jack-knife procedure with bin size of 200 iterations. The exceptional configurations are 
discarded with the same criterion as defined in Eq. (|15|l . Results of hadron masses are collected in Tables lXXXVI - 
IXLI in Appendix A. 



B. Static quark potential 



We calculate the static quark potential in order to determine the Sommer scale p9| which we use in our 
analysis of hadron masses. For this purpose, the temporal Wilson loops W{r, t) up to t = 16 and r = 
are measured both in full and quenched QCD simulations. Wc apply the smearing procedure of Ref. |36l | up to 
twelve steps and the measurements are carried out every four steps. 

The static quark potential V{r) is determined from the correlated fit of the form 

Wir,t) = C(r) exp [-F(r)t] . (17) 

We take the fit range [iminj^max] — [3,7] in all simulations in full and quenched QCD by inspecting the t 
dependence of the effective potential 

Feff (r, t) = In [W{r, t)/W{r, t + 1)] . (18) 

Examples of VcS are plotted in Fig. For each r, the number of smearing steps is fixed to its optimum value 
at which the overlap to the ground state C(r) takes the largest value. 

As shown in Fig[7| we do not observe any clear indication of the string breaking. Therefore V{r) is fitted to 
the form 

V{r) = Vo-- + ar + SV{r), (19) 
r 

where 5V{r) is the lattice correction to the Coulomb term calculated perturbatively from one lattice-gluon 
exchange diagram |37| 

SV{r) = -5, (cir) - , (20) 
Gir) ^ 4.V 3^"^[^J] . (21) 



The Sommer scale tq defined through |2E 



dV{r) 



'° dr 



1.65 (22) 



is then determined from the parametrization of the corrected potential V{r) — SV(r): 



/ 1.65 ~ a 

.„ = V^^- (23) 

The lower cut of the fit range in Eq. H19|l . rmin, is determined by inspecting the r^i^^ dependence of rg. We 
observe that is relatively stable for Vn^in & [V2, 2^2] as shown in Fig|S| With r,„in < V2, x^/dof takes 
an unacceptably large value due to the violation of rotational symmetry, while a becomes ill-determined with 
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''mill > 2-\/2- We therefore take rmin = -s/S- While the Tmax dependence of tq is rather mild, the covariance matrix 
becomes ill-determined with rmax greater than 9-^2 on 20'^ x 48, 7-^2 on IG'^ x 48 and 6\/2 on 12^ x 48. We 
therefore fix rmax to these values. 

We repeat the fits, Eqs. (|17|) and (|19|) . with other choices of the range: ^min = 4 or r^ain G [V2, 2\/2]. The 
largest deviations in the fit parameters and tq are included into their systematic errors. Other systematic errors 
due to the choice of imax, the optimum number of the smearing steps and Tmax are small and ignored. Fit 
parameters in Eq. (|19|l and tq are summarized in Tabic IITI for full QCD, and in Tabic IITTI for quenched QCD. 



C. Autocorrelation 



The autocorrelation in our full QCD data is studied by calculating the cumulative autocorrelation time 

At„ 

-I- 

2 



tS^-(AW) = J + E /'o(At), (24) 
where po{t) is the autocorrelation function 



ro(At) = ((0(i) - (O)) (0(i + Ai) - (O))) (26) 



and we take Aimax = 200. 



In Table HVl we summarize rg^™ in the A-HMC and S-HMC simulations on the 20^ x 48 lattices for three 
quantities: i) the plaquette which is measured at every trajectory, ii) the PS meson propagator at t=12, iii) the 
temporal Wilson loop with (r, t) = (5, 4). The results do not show any systematic differences in r^"™ between the 
A-HMC and S-HMC runs. The plaquette shows the largest autocorrelation with '''pi™ = 10-30, which is similar 
to those found in the UKQCD simulation using the same lattice action and similar simulation parameters. 
We obtain smaller values of r^"" for the other two quantities. This is contrary to a naive expectation that 
these long-distance observables have a longer autocorrelation than the local quantity like the plaquette. This 
suggests that the size of noise arising from short correlation modes is larger than that of the longest mode in 
these observables and our statistics are not sufficient to extract t^"'° of the longest but weak mode. 



The statistical error including the effect of autocorrelation is given by ^2t™™ times the naive error. There- 
fore, the above observation tells us that the bin size in the jack-knife procedure of 60 HMC trajectories or larger 
is a safe choice to take account of the autocorrelation in our data. 

The bin size dependence of the jack-knife error of hadron masses and the static potential is plotted in Figs. El 
and 1101 We use errors obtained from uncorrelated fits because, with large bin sizes, the number of bins would 
not be sufhciently large to determine the covariance matrix reliably. For both hadron masses and the static 
potential, the jack-knife error reaches its plateau at bin size of 50-100 trajectories, which is roughly consistent 
with the above estimate from r^"™. The situation is similar on smaller volumes 16'^ x 48 and 12^ x 48. We 
therefore take the bin size of 100 trajectories in the error analysis in full QCD. 

We also investigate the bin size dependence of the jack-knife error in quenched QCD. As shown in Figs. 1111 
and 1 121 the bin size of 200 iterations is reasonable. 

Another point of interest is the sea quark mass dependence of the autocorrelation. A natural expectation 
is that smaller sea quark mass leads to a larger correlation length and hence a longer autocorrelation. This 
expectation is supported by the CP-PACS observation in Ref. 11], where they used the RG-improved gauge and 
clover quark actions. However, our result of Tpj"™ in the S-HMC simulations, which is determined more precisely 
than that for the A-HMC due to the higher statistics, shows the contrary sea quark mass dependence: Tpj™ 
decreases as the sea quark mass decreases. This is consistent with the UKQCD's observation in Ref. ^9]. We 
also note that Tpj"™ in our simulations is much larger than in the CP-PACS 's runs particularly at the heaviest 
sea quark masses mps,sca/'7iv.sca — 0.8 (i^sca^ 0.1340). 

In our determination of non-perturbative csw at /3 = 5.2 (23|, we find that the expectation value of the 
plaquette varies rapidly around Kgca. — 0.132, where the plaquette shows the strongest autocorrelation in the 
investigated region K E [0.100,0.136]. Since such a behavior, somewhat similar to a phase transition, is not 
observed at higher /3, we consider the unexpected behavior of the plaquette to be an artifact due to finite lattice 
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spacing. This artifact is probably absent or well suppressed with the CP-PACS's choice of the improved actions. 
At sufficiently small lattice spacings, we then expect that Tpj"™ shows the natural sea quark mass dependence, 
namely larger Tpj"™ for lighter sea quark masses. We also expect that, even at 13 — 5.2, Tpj"™ will increase if the 
sea quark mass becomes sufficiently small compared to that corresponding to .ftTsoa^O.f 32. 



IV. FINITE-SIZE EFFECTS 



Finite-size effects (FSE) are one of the major sources of systematic errors in lattice calculations. Since our 
largest volume size ~ (1.8 fm)'^ is still not so large, it is important to check FSE in our data. We discuss how 
much FSE is present in our data on the largest lattice using data on three spatial volumes 12^ x 48, 16^ x 48 
and 20^ x 48. 

In Figs. I13 H15I we plot diagonal data of hadron masses as a function of the spatial volume inverse. For 
'Tips,soa/'7iv,soa^0.7, including the quenched case, hadron masses obtained on the 16'^ x 48 and 20'^ x 48 lattices 
are consistent with each other within two standard deviations. On the other hand, hadron masses decrease 
monotonously up to V = 20'^ at the lightest sea quark mass corresponding to ?Tips,soa/"^v.sca — 0.6. 

The magnitude of FSE also depends on the valence quark mass. Figure El shows the valence quark mass 
dependence of the relative mass shift between the two larger lattices for PS mesons and octet baryons. We 
observe that, except at the lightest sea quark mass, the mass shift is at most a few percent level in the whole 
range of the simulated valence quark mass. The situation is similar for vector mesons and decuplet baryons. 
Therefore, we conclude that the size of FSE on our largest lattice is small over our range of valence quark masses 
down to the second lightest sea quark mass. 

The mass shift is non-negligible at the lightest sea quark mass. While the magnitude is of the order of a few 
percent for the heaviest valence quarks, it clearly increases as the valence quark mass decreases. 

We consider that the observed FSE is caused by valence quarks wrapping around the lattice in spatial 
directions (namely squeezing of hadrons into the small box) rather than wrapping of virtual pions. As shown in 
Fig. El the magnitude of FSE caused by the effects of virtual pions (long dashed line) [sll given by mps(i) — 
mps{L — oo) ^ exp [—mps{L — oo)L] with L — N^a is too small compared to observed effects. 

A qualitative understanding of the observed FSE is as follows. The wrapping of valence quarks is suppressed 
by the center Z{i) symmetry in quenched QCD p^ . In full QCD, Z(3) symmetry is broken by the wrapping 
of sea quarks in the spatial directions, whose magnitude increases toward lighter sea quark. A possible reason 
why FSE is significant only at our lightest sea quark mass would be that the Z(3) breaking turns on rather 
quickly around the lightest sea quarks. 

The enhancement of FSE toward the lighter valence quarks leads to a decrease of the slope dmhad/drripg in 
Fig. 1171 and, hence, underestimation of hadron mass splittings, such as the K*-K hyperfine splitting. The mass 
splittings are expected to be increased by sea quark effects, since these are underestimated in quenched QCD 
as well established in Ref. 4]. Therefore, FSE makes sea quark effects less clear. It is crucial to check how large 
FSE is in our hadron mass data at the lightest sea quark mass on the largest lattice. 

Figures [Ql and 174] show that the volume dependence of our data is well described by a power law 

'71had(^) = mhad(i = oo) + c/L^ (27) 

as found in Ref. using the KS fermion. The relative size of FSE on the largest lattice 

TO(L==20a)-mhad(i = oo) 

Am = ^ 28 

m(L = 20a) ^ ' 

is estimated from this ansatz and is plotted in Fig. El We find that, for PS and vector mesons. Am is about 
5 % for diagonal data and is reduced to a few percent at _R'vai = 0.1350, which roughly corresponds to the strange 
quark mass m^. It is expected, however, that the volume dependence (|27|l turns into a milder form exp[— mpgL] 
for sufficiently large volumes. The actual size of FSE should be smaller than the above estimation, say, a few 
percent. Since this is smaller than the typical size of quenching errors, which is 5-10%, we consider that the 
examination of sea quark effects is feasible in the meson sector, particularly in strange meson masses. 

Finite-size effects are more pronounced for baryon masses as observed in Fig. 1181 For diagonal data. Am is 
roughly comparable with typical quenching errors in the baryon spectrum of the order of 5-10 %. Sea quark 
effects in the light baryon masses, such as m^r and m\, may become unclear by the contamination of FSE. 
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We note, however, that Am decreases for heavier valence quark masses. The examination of sea quark effects 
becomes more feasible for strange baryon masses like toh and niQ. 

Figure [THI shows Am for the diagonal data of the octet baryon mass as a function of Ng. We find that the size 
of A^s~30, which corresponds to LKi2.7 fm, is required to suppress Am to a few percent level. The required 
size becomes slightly smaller, Lk. 2.4 fm, for the valence quark mass around rUs- These sizes are larger than 
our largest spatial lattice size L~1.8 fm. Further simulations on such large lattices will be needed to obtain a 
definite conclusion on sea quark effects in the baryon spectrum. 

Tables UTI and imi show that FSE in tq are much smaller than in hadron masses. While the central value on 
20^^ is systematically higher than that on \& in both full and quenched QCD, the difference is about 1 % and 
not significant with the accuracy of our data. The size of FSE is small, namely a few percent level, even on 12^. 
Therefore FSE in on the largest lattice can be safely neglected both in full and quenched QCD. 



V. CHIRAL EXTRAPOLATION 



The hadron spectrum in full QCD is calculated using hadron masses measured on the 20'^ x 48 lattice. This 
requires a parametrization of the mass data as a function of sea and valence quark masses in order to extrapolate 
(up-down) or interpolate (strange) quark masses to their physical values. We make this parametrization by 
combined fits to masses of a given hadron at all sea quark masses. We test the following two methods for the 
combined fit: 

A The effective lattice spacing, determined from for instance, may vary as a function of the sea quark 
mass. In order to separate this effect from the physical quark mass dependence, we carry out the 
chiral extrapolation using dimensionless quantities such as (ro(ifsoa)'Tihad(^sea; -f^^vai,!, -f'^vai,2)), where 
'Tihadl^sca; -f^vai,!, ^vai,2) represents the measured hadron mass composed of valence quark masses cor- 
responding to -R'vai 1 and Xvai 2 on the gauge configurations generated at i^soa- We refer to this way as 
"method- A". 

B It is also possible to fit hadron masses in lattice units, as was done by the SESAM and the CP-PACS [Tl| 
Collaborations. We call this "method-B". 

A detailed description of the two methods will be given in Sees. IV Al and IV Bl They should yield a consistent 
hadron spectrum, since fit forms in method-B can be reproduced from those in method-A by expanding rg as 
a function of sea quark mass. This consistency is examined in Sec. IV Bl 



A. chiral extrapolation using ro (method-A) 



Chiral perturbation theory (ChPT) "4^ provides a guide to obtain a controlled chiral limit of hadron mass 
data. For the quark mass dependence of diagonal data of the PS meson mass, ChPT predicts the presence of 
logarithmic singularities. At the one-loop level, the ChPT prediction reads 



'PS 



l + ^v\n[v]+Av (29) 



2Bom, Nf 

where y = 2BQmq / {4:Tt f)"^ and A is a linear combination of the low energy constants of ChPT Lagrangian 
: A = {2as — cx^) + Nf {2ae — 0:4). The mass ratio on the left hand-side, mpg/mq, is plotted as a function of 
niq (X y in Fig. 1201 As we already reported [21I |23| , our data show no hint of the curvature predicted by the 
chiral logarithm. The fit of Eq. H29|l . assuming / to be a free parameter, gives / ~ 6 GeV, which is much larger 
than its experimental value 93 MeV. On the other hand, the fit gives an unacceptably large x^/dof = O(IOO), 
if we fix / to the experimental value. 

A similar test using formulae from partially quenched ChPT (PQChPT) "4T| also shows that the coefficient of 
the chiral logarithm term obtained from our data is much smaller than the prediction from PQChPT [2lll2^l4^ . 
A possible reason for the absence of the chiral logarithm is that the sea quark mass in our simulations is still 
too large and higher order corrections of ChPT should be included to describe the data. 



In this study, therefore, we use simple polynomial fitting forms in terms of the quark mass for the chiral 
extrapolation. The systematic error due to the chiral extrapolation is estimated by testing several different 
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polynomial forms. However, the inconsistency between our data and ChPT suggests that the extrapolation 
may have larger uncertainty than this estimation. This point will be examined in detail in a separate paper 

m 

Since the sea quark mass in our simulations is not so small as discussed above, it is important to check the 
convergence property of the polynomial expansion of hadron masses in our range of the sea quark mass. We 
carry out both quadratic and cubic chiral fits to diagonal data of PS and vector meson masses: 

^Ps(-^sca7 -^scai -^sca)) ^ -^diag /^q,sca H~ ^diag Mq,sca ^~ -^diag Mq,sca' 

(30) 

ro{Ksca.) mY{Ksca.] Ksca, Ksca) — ^diag + ^diag MPS.sca + C^ag MpS,soa + -^diag MpS.soa' (31) 

where ^q,diag and /ips.sca are the quark mass defined through the vector Ward identity (VWI) and the PS meson 
mass normalized by tq: 

Mq.sca — ^0 (-^sca) ^q.seai (3^) 
1/1 1 



MPS, sea — ,sca r, (34) 

Wps,soa = mps{Ksca;Ksca.,Ksca)- (35) 

Fit parameters and x^/dof are collected in Tables IVl and IVTI The coefficient of the cubic term is small and 
consistent with zero for both Eqs. (|30() and (|31|l . Consequently, the quadratic and cubic fits show a good 
consistency with each other in the whole range of the quark mass and toward the chiral limit, as seen in Fig. 1211 
These observations suggest that the polynomial expansion up to the quadratic order is sufhcient to describe the 
quark mass dependence of our data of meson masses in the method-A. 

We carry out a combined fit to PS meson masses as a function of sea and valence quark masses using the 
quadratic form 

— rPS ,, _i_ rPS ,, I r'PS,,2 I /-"PS ,,2 I /-"PS,, ,, fQC\ 

— -Ds Mq,soa i- -Dv Pq,val + L-g Mq,sca + Pq,val + ^sv A'q.sca A'q.val , I^ODj 

where /iq,soa is defined in Eq. H32|l and 

Mq,val = ro{Ksca.) rriq^ys.i, (37) 
Wq,val = ^ ('nq,val,l + TOq,vaL2) , (38) 

^q.aM = ^^-^) (-1,2). (39) 

The presence of the monomial term in mq_soa means that the PS meson mass does not vanish in the chiral limit 
'Tiq.vai = for non-zero values of rriq^soa- This is because the value of iiTvai where the PS meson mass vanishes 
depends on the sea quark mass due to explicit violation of chiral symmetry with the Wilson-type quark action. 

We employ uncorrelated fits in the combined chiral extrapolations although the data with the same sea quark 
mass are expected to be correlated. Therefore, the obtained x^/dof can be considered only as a guide to judge 
the quality of the fit. Figure [23 shows that this fit form describes our data well. Parameters of the fit are 
summarized in Table IVIII We note that Kc determined from the diagonal fit Eq. (|30|l and the combined fit 
Eq. H36|l are consistent with each other, as they should be. 

The most general quadratic fit ansatz for the PS meson masses should include an additional cross term 

{roiKsca) TOPS (Xsca;ifvaU,^vaL2))' = "r.h.S of Eq. (EH)" 

/^q,val,l Mq,val,2; 

(40) 



where /iq,vai,i = '"o(^sca) TO,q^vai,i- Howcvcr, the coefficient is small as shown in fit parameters in Table IVIII and 
hence does not change the hadron spectrum. We use this fit to estimate the systematic error due to the choice 
of the fitting function. 

For the vector meson, wc find that the following form describes our data well 

roiKsca,) mv(i^soa; ^val.l, ^i^val,2) = A^^ + B'^ MPS.soa + MPS.val + MPS.sea MPS.val, (41) 
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where /ips.soa is defined in Eq. H34|) and 

MPS.val = (7'o(J'^sca)"lps,val)^, (42) 

"^i■S,val = ^ (™PS,val4 + "^PS,val,2) : (43) 

Wps.val.i = Wps(i^soa;-f'^val,i,^val,i)- (44) 

For a more general fit of form 

raiKsca) my{K,^a] K^aiA^K^aia) = "r-^-s of Eq. gTJ" + C^Mpg.eca + C^Mls.vaP (45) 

the additional parameters and are not well-determined as seen in Table fVlTTl We use the former fit, 
which is shown in Fig. to calculate the hadron spectrum, and the latter to estimate systematic error of the 
chiral extrapolation. 

We also carry out a partially quenched fit to vector meson masses at each sea quark mass. We use a linear 
form, which is obtained from Eq. |(1T1) by dropping all terms describing the sea quark mass dependence 

mv(i^sca; ^vaLli ^val,2) = ^PQ + -B^Q "lpS,val- (46) 

Parameters given in Table Hxl are used to calculate the J parameter at each sea quark mass. 



The chiral extrapolation of octet baryon masses is carried out using a quadratic form based on the leading 
order prediction of ChPT 43], which was also used in Ref. We carry out the simultaneous fit to the S- 

and A-like octet baryon masses using the functions 

^ (-^sca ) moct,T.{Kaca.] -f^val , 1 , -f^val , 2 , ^val,2) 
= A° + i?f ^PS.soa + {F^ - F>^) AiPS.vala + 2F°^ps,val,2 + C's°A'PS,sca + (C^ + C^?'^) ^4>S,val^ 

Mps ,val,2 MPS, sea MPS,val,l MPS, sea MPS,val,2, (47) 

''o(^SGa) '71oct,A(^soa; -f^val , 1 , -?^val , 2 , -K'vaL2) 

= A'^ + B^flps.sca + (^F^ H MPS,val,l + 2 ^F^ — 3^°) /^PS,val,2 

Mps ,sca ,val, 1 MPS,val,2 + (^sv + ^sv' ) MPS,sea MPS,val,l 

"^(^sv ^sv ) MPS, sea AiPS,val,2, (48) 
where /ips^vahi = (?'o(-^sca) mps,val,i)^- 

The decuplet baryon masses are well described by the following form 

roiKsca) m^cciKsca] Kva\,lT Kva\,2, Kva\,2) ^ + i?f A'PS.sca + MPS.val + C*^ /^PS.sea A*PS,val, (49) 

where /xps,vai stands for the average of three valence quark masses 

AiPS,val = {rQ{Ksca)mps,-va\f , (50) 
"^PS,val = 2 (™PS,vaLl + 2mpg ^^1,2) ■ (51) 

Figures El and 123 show the fit for octet and decuplet baryon masses. Parameters are summarized in Tables Ixl 
and IXII We also test the following forms to estimate systematic error of the baryon spectrum due to the choice 
of the fitting form 

rQ{Ksca) rnoct{Ksca\ K^a\S, Kval,2, Kval,2) 

= "r.h.s of Eqs. 63) and + C° MPS,vai,i /^PS,vai,2, (52) 

raiKsca) ™dcc(-ft^sca; -^^val , 1 , ^val , 2 , ^val,2) 

= "r.h.s of Eq. 63" + C^^^l ^ _ + C^fils.^.^ (53) 

We carry out the chiral extrapolation of rQ{Ksca) in order to determine vq at the physical sea quark mass, 
which is required to calculate the hadron spectrum in the method-A. We use a linear form 

1^A,, + |^. (54) 
As seen in Fig. 1261 and Table Ixlll this fit describes our data well and gives a reasonable value of x^/dof~1.5. 
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B. chiral extrapolation in lattice units (method-B) 



In order to study the convergence properties of polynomial fit forms in the method-B, we carry out quadratic 
and cubic diagonal fits to PS and vector meson masses 

"^PS (^soa; -K^sca, Ksc&Y — ^dTag "^q,sca + C'di'ag "^q.soa + ^cUag "^q.soa) (55) 
mv{Ksca; Ksca.: Ksca.) — ^dYag + -^dYag "^PS,sca + ^'dTag 'Tips, sea + -^dYag '^PS.sea- (56) 

Fit parameters are summarized in Tables IXIIII and [XIVI 

Fit curves of the quadratic and cubic fits to vector meson masses are shown in Fig. 1271 While quadratic and 
cubic fits describe our data reasonably well at quark masses used in the simulation, they develop a deviation 
toward the chiral limit and for heavy quarks. 

Figure 1281 compares the relative magnitude of the linear and quadratic terms in the quadratic diagonal fit 
to PS meson mass in the method-A and B. As the quark mass increases, the magnitude of the quadratic 
contribution in the method-B increases more rapidly than in the method-A; it is no longer a small correction at 
the simulated quark masses (mq^diag — 0.015-0.055). A similar situation is observed in the chiral fit to vector 
meson masses as shown in Fig. |2ni 

We come to conclude that the chiral expansions of meson masses in lattice units, Eqs. H55|) and (|55|) . have 
poor convergence properties compared to those in unit of tq, Eqs. H3U|I and (|31|l in the method-A, and the cubic 
term should not be ignored in the method-B. We directly confirm this point in Fig. 1301 where the fit results 
for vector mesons from the method-A are converted to lattice units using Eq. H54I) and compared with the fits 
of method-B. The cubic fit in the mcthod-B shows a good consistency with the quadratic fit in the method-A, 
while the quadratic fit in mcthod-B does not. 

The combined chiral fit including cubic terms is not very stable because it contains a number of free param- 
eters. In this study, therefore, we do not use the method-B to extract the physical hadron spectrum. 

Before we turn to details of the determination of the hadron spectrum, let us make additional comments on 
the failure of the method-B with our data. Figure 1^ shows the chiral fit of tq as a function of Wpg 

I'oiKsca) — + ^ro "^PS.sca + ^ro "^PS,sca + ^ro '^PS,soa- (57) 

Fit parameters in Table IXVl show that a large contribution of higher order terms is present also in this fit. By 
substituting this parametrization of tq to Eqs. H3U|) and Ij^lfl (diagonal fits in the method-A), the large higher 
order corrections appear in Eqs. (|55() and (|56|l (diagonal fits in the method-B). Conversely, why the method-A 
works well is that large contributions of higher order terms in hadron masses and tq cancel with each other at 
least partially. 

We note that the method-B works well in the CP-PACS's study where they took similar simulation 
parameters but with different lattice actions, namely the RG-improved gauge action and the tadpole improved 
clover quark action. We compare the CP-PACS data of romy at (3 = 2.2 and ours in Fig.|211 A good consistency 
in the whole range of the quark mass suggests that two groups' data are in the scaling region. However, we 
find that the CP-PACS data of ro in lattice units show much milder dependence on the sea quark mass than 
ours. This is the reason why the method-B works well in the CP-PACS study, but does not with our data. It 
is, of course, not surprising that different lattice actions lead to different sea quark mass dependences of hadron 
masses and tq in the lattice units. However, as discussed above, the much stronger dependence with our choice 
of the lattice action is practically problematic, if one carries out the chiral extrapolation in lattice units. 



C. calculation of hadron spectrum 



The hadron spectrum at the physical quark mass is determined as follows. The pion and p meson masses 
normalized by tq are determined by tuning their ratio (roTO^)/(ror7ip) to its experimental value, i.e., by solving 
the equation, 

(roTOTr) ''^7r,exp /KQ\ 

, [0?i) 



AV + (BY + BY) (rom,)2 + Cl{r^m^Y 
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where we denote the experimental value of hadron mass by mhad,exp, and rg represents the Sommer scale at the 
physical sea quark mass. The hopping parameter corresponding to the physical light quark mass, Kud^ is fixed 
by solving 

{ r.h.s of Eq. with ii:sea = Ji^vai,i = -^vai,2 = ^S^«d} = (rom^)^. (59) 
Then we determine rg, which is required to convert r^rriTj and r^mp to TOtt and m^, from Kud and Eq. H54|) . 
We test two meson mass inputs to fix the strange quark mass: 



• In the first method, we use the kaon mass as input. The hopping parameter corresponding to the strange 
quark mass, K^, is determined by solving 



r.h.s of Eq. _ m,K,c^p ^^^-^ 
romp mp^c^p ' 

where we set Ksca = ifvai,! = Kud and Js'vaL2 = Kg in the r.h.s of Eq. Then the mass of the 

meson, that is an unphysical ss PS meson, is determined from Eq. H36II . and used to calculate strange 
vector meson and strange baryon masses. We refer to this meson mass input as if-input. 

• In the second method, we use the 4> meson mass as input assuming that it is a pure ss vector meson. The 
?7s meson mass is fixed from 

r.h.s of Eq. jH)) _ m^.cxp ^^^^ 



romp mp^cxp 

where we set mps^scn — m-^ and TOps,vai,i = '7ips.vai,2 = "^J73- We determine Kg from to^^ and Eq. H36|l . This 
input is called </>- input. 



The full spectrum of non-strange and strange hadrons is determined by substituting K^d, m„, Kg, m^^ and 
ro to Eqs. 1)36(1 . ((41|l . 147(1 . ((48(1 and 149(1 . We use the lattice spacing a determined from mp to convert the 
hadron masses in lattice units to those in physical units. We note that this estimate of the scale a is subject to 
a systematic uncertainty due to the use of the polynomial fitting forms for the chiral extrapolation. However, 
if we use ro as the input to set the scale, we obtain a consistent result for a within errors. The results of Kud, 
Kg and are collected in Table IXVH 

We repeat the above analysis using each of the alternative fit forms Eqs. ((4U(I . ((45(1 ((52(1 and 1(53(1 . The largest 
deviation in the hadron spectrum among these analyses is taken as the systematic error due to the choice of 
the chiral fit forms. 

For the chiral extrapolation of ro, we find that an alternative form 

ro{K,,.)=Ko+B'rAs,s.. (62) 

also describes our data well. However, the hadron spectrum calculated using this fit is completely consistent 
with those using Eq. 1(54(1 . We therefore ignore the systematic error due to the choice of the fit form Eq. ((54(1 . 

The systematic error of the measured value of ro{Ksca) leads to an additional uncertainty in the result of the 
hadron spectrum. We perform the calculation of the spectrum with ro{Ksca) shifted by its systematic error at 
one value of i^soa- This calculation is repeated for all Ksea. and the largest deviation in the spectrum is included 
into the systematic error. 



D. chiral extrapolation in quenched QCD 



The chiral extrapolations in quenched QCD are performed using fit forms which are obtained from those used 
in the full QCD analysis by dropping all terms describing the sea quark mass dependence. Namely, fitting forms 
for meson masses are 

mps(i^val,l,-ft^val,2)^ ^ TOq^val + ^ ^^^^j, (63) 

mv(i^val,l,^i^val,2) = + ,,^1. (64) 
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The following fitting forms are used for baryon masses 

^PS,val,2' 



+ (CO-C'°^^)m4g (65) 




moct.A(ifvaU,^val,2,ifval,2) = ^q + | ^q"" + ^ ) ^^^..1,1 + 2 (^^q" - m|s,val,2 

+ + C-O'^) ml,s,^,,^, + (C-o - C"^^) m4s_^^,_2, (66) 

nT-doc(^val,l>^val,2,^val,2) = + 5° mlg^^^i . (67) 

These forms fit to our data very well as shown in Fig. 133 Fit parameters are summarized in Table IX VIII The 
hadron spectrum is calculated in an analogous way to that for full QCD. Resulting values of Kud, Kg and 
are summarized in Table IXVlIII 



VI. SEA QUARK EFFECTS ON HADRON SPECTRUM 

A. Sea quark effects at simulated quark mass 

Figure |M1 compares the valence quark mass dependence of the vector meson mass at each sea quark mass 
in full QCD and in quenched QCD. We observe that quenched data have a significantly smaller slope than 
experimental points. This leads to the underestimation of the K*-K hyperfine splitting in quenched QCD. 

The slopes in full QCD data are clearly larger than in quenched QCD, and increase for decreasing sea quark 
mass. This is reflected in a negative value of C^^ for Eq. H41|l in Table IVIIII This sea quark effect leads to a 
better agreement of the meson spectrum in full QCD with experiment than in quenched QCD. 



The J parameter [4J| defined by 



J — my 



dnipg 



(68) 

mps/mv=mj<:/"ijf • 



is useful to quantify the sea quark effect. Numerical results of J calculated from the partially quenched chiral 
fit Eq. H46|l are given in Table IXIXI In Fig. we plot J in full QCD (filled circles) as a function of the sea 
quark mass together with the quenched result (open circle in the right panel). We observe that J in full QCD 
is close to the quenched value at heavy sea quark masses corresponding to mps/my > 0.75 and increases as the 
sea quark masses decreases. 

In the same figure, we also plot J reproduced from the combined chiral fit Eq. H41|) (dashed lines). The 
result is consistent with J from the partially quenched fit, as it should be, and shows a similar sea quark mass 
dependence. We observe that J extrapolated to the physical sea quark mass is closer to the phenomcnological 
value 13 

J ^ niK^ ~ = 0.48(2) (69) 
mj^ — raj^ 

than in quenched QCD. 

Figure 1331 also shows J calculated from our results of the meson spectrum (see Table IXX|) using the above 
alternative definition Eq. I|69() (filled square). The result is in good agreement with other determinations, 
showing the magnitude of the sea quark effect in J to be stable against the definition of J. 



In Sec. IIVI we pointed out that FSE decreases the slope dmy / drripg. This is confirmed numerically in 
dmv /drripg determined from the partially quenched chiral fit Eq. at the lightest sea quark mass: dmy /drripg 
is 0.906(14) on 20^, 0.814(88) on 16^ and 0.68(30) on 12^. The slope would be larger if we increase the spatial 
size beyond 20''. Therefore, the observed sea quark mass dependence of the slope and J is a genuine effect of 
dynamical quarks and not an artifact of FSE. 

A similar effect of sea quarks can be found in decuplet baryon masses as shown in Fig. 1361 However, a 
significant deviation still exists in the slope between full QCD data and the experimental spectrum. We consider 
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that a larger slope in full QCD is still partly masked by FSE on the 20'^ volume; volume as large as 30^ would 
be needed to reduce FSE to a few % level as discussed in Sec. IIVI 

We emphasize that the evidence of sea quark effects observed in this subsection does not suffer from possibly 
large systematic errors due to the chiral extrapolation and the choice of inputs to fix the scale and quark masses: 
the increase of the slopes, dmy / drripg and c?TOdcc/'7ipg, is observed without any chiral extrapolation and inputs. 
The sea quark mass dependence of J is obtained by a short extrapolation or an interpolation to a relatively 
heavy valence quark mass corresponding to mps/mY — rriK /friK* ■ 



B. hadron spectrum 

The meson spectrum in full and quenched QCD is summarized in Table IXXl Since our fitting functions to 
vector meson masses, Eqs. ()41|l and (|64f) . are linear in terms of the valence quark mass, m^* with (ji- input 
equals (rrip^exp + '^i^.cxp) /2 in both full and quenched QCD. The deviation of this value from the experimental 
mass mif.^exp is only 0.2%. 

A clear difference between full and quenched QCD is observed in other meson masses as shown in Fig. IS7I 
While the quenched meson spectrum shows a significant deviation from experiment, sea quark effects reduce 
the deviation by about 40%. This closer agreement of the meson spectrum with experiment is a consequence 
of the sea quark effects observed in the previous subsection. 

In Fig. I2HI mK* with the if-input in full QCD is compared with the CP-PACS results obtained with the 
RG-improved gauge and clover quark actions |ll|. We observe that our rriK' is consistent with the CP-PACS' 
result at a similar lattice spacing, and is at the lower edge of their estimate in the continuum limit. In the same 
figure, we also make a comparison in quenched QCD with the CP-PACS results obtained with the plaquette 
gauge and the Wilson quark actions in Ref. Two groups' results show a good agreement with each other. 
These observations suggest that the scaling violation is small in our data both in full and quenched QCD and, 
hence, the closer agreement of the meson spectrum in full QCD with experiment is a genuine effect of sea quarks. 

The baryon masses in full and quenched QCD are listed in Table IXXII These masses are compared with 
experiment in Fig. 1391 For heavier baryons, such as S, S and 17, full QCD results show a closer agreement 
with experiment than in quenched QCD. The sea quark effect is, however, less clear for lighter baryons. This 
is partly due to FSE in full QCD data which is more pronounced for lighter valence quarks. 

In Fig.Hni and with the /-C-input are compared with the CP-PACS results j3,[l3- While the fuh QCD 
results of two groups show a reasonable agreement with each other, the CP-PACS results in the continuum limit 
in quenched QCD are systematically smaller than ours. This suggests that our quenched data has non-negligible 
scaling violation, which is another source making sea quark effects less clear. Therefore, further investigations 
of FSE in full QCD and scaling violation in quenched QCD are required to obtain a clear conclusion on sea 
quark effects in the baryon spectrum. 

We now turn to theoretical predictions which can be derived from our data. The first is the mass of the 
rjs meson, for which our full QCD data predict to^, = 0.6948(3) (4-8/— 1) (+2) GeV with the iiT-input, and 
0.7381(46)(-t-57)(-|-40/— 46) GeV with the 0-input, where the first error is statistical, and the second and third 
ones are due to the choice of the fitting form and a systematic uncertainty of tq. These results are to be 
compared with those in quenched QCD, 0.6988(9) GeV (X-input) and 0.7719(58) GeV (0-input). While the 
values themselves do not differ by going from quenched to full QCD, the difference between the two inputs is 
reduced by about 40% in full QCD. This reflects the closer agreement of the meson spectrum in full QCD. 

Another interesting prediction is the physical value of rp. Our full QCD simulation gives — 
0.497(6)(~9)(-|-ll/— 12) fm, where the meaning of three errors are the same to those of m,,^. We note that 
this is close to the phenomenological estimate in the original paper |29|. ro=0.49 fm. The quenched simulation 
gives 0.5702(75) (50) fm, where the first error is statistical and the second comes from the systematic uncertainty 
of the measurement. About a 14 % difference between full and quenched QCD arises from the following two sea 
quark effects. One is the difference of the physical value of tq itself due to the change of the shape of the static 
quark potential. The other source is the reduction of the quenching error in nip in full QCD, which is used to 
fix the lattice scale. 
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VII. DECAY CONSTANTS 



The PS meson decay constants are calculated using the fourth component of the improved axial vector current, 
which is defined by 

AT = Ai + CA^iP (70) 

with the symmetric lattice derivative A4. We extract the amplitude of the (A™^(t)P(0)'l') correlator, C^^, by 
the correlated fit of the form 

(AfP(OP(O)t) = {exp [-mps t] - exp [-mpg {N^ - t)]} (71) 

with mps fixed to the results given in Appendix A. We use the local operator for A™^, while the double 
smearing is applied to P. The amplitude of the (P(t)P(O)^) operator with the doubly smeared source and the 
sink operators, Cp^, is extracted assuming a single hyperbolic cosine form. 

The renormalized decay constant is calculated by 



/ps = 2ifZ^(l + &Amq,vai)C^%/ (72) 

y mpsC^'' 

Since non-perturbatively determined values for Z4, 64 and C4 are not available for two-flavor QCD, we adopt 
one-loop perturbative values in Refs. [iR Wn. 14^ lifll with the tadpole improvement. We calculate 
ap(3.40/a) from the plaquette average {Up) according to [SU Is^ 

47r 

-In [(C/p)] = —ap(3.40/a){l- (1.1905 - 0.22667V/) ap}. (73) 

Then, Qp(3.40/q) is evolved to the optimum scale (q*^^^ =1.803/a for Z^, ql^ =2.289/a for 6^ and q*^ =2.653/a 
for CA [53) using the universal two- loop beta function and is used as the expansion parameter of tadpole 
improved perturbation theory. 

The consistent chiral extrapolation of the decay constant should include the chiral logarithmic term as pre- 
dicted by ChPT 40]. However, our data do not show the characteristic curvature of the chiral logarithm as 
discussed in Sec IV Al fand also in Refs. 0,11^). We therefore use the following polynomial form for the chiral 
extrapolation leaving the problem of the chiral logarithm and associated uncertainty for future publication |42| 

?'o(-f^sca) /ps(-f^sca; K , K ^2) ~ A^ + BI. /ips.sca + B\, /ips,val + C'f^ AiPS.sca MPS.val- (74) 

The fit is plotted in Fig.^Jwith parameters summarized in Table IXXIII Pion and kaon decay constants, /^r and 
fx, are calculated by tuning /ips.sca to {r^niT^)'^ and substituting (rom^)^ or {roruri^)'^ for //ps.vai.i (* = li2) in 
MPS.vai- The systematic errors due to the choice of the fitting function and the uncertainty of rg are estimated 
in a way similar to those described in Sec. IV CI In the estimation of the former error, we use 

ro(i^sca) /Ps(i^sca; i^vaU , i^val , 2 ) = "r.h.S of Eq. d" + Cf/i|s.sca + MpS.val (75) 

and Eqs. H40I) and (|45|l as the alternative fitting functions for the chiral extrapolation. 
In quenched QCD, we use the chiral extrapolation form 

/Ps(ifval,l,^val,2) = < -f m|s^val • (76) 

and obtain parameters summarized in Table IXXIIII For Za, bA and ca, we test the one- loop perturbative value 
and the non-perturbative one in Ref . [s^ . 

Our results of the decay constants are summarized in Table IXXIVI A comparison between full and quenched 
QCD results obtained by the one-loop matching is made in Fig. 021 We observe that /tt is consistent with the 
experimental value within two standard deviations in both full and quenched QCD. While fx in quenched QCD 
is significantly smaller than the experimental value, the deviation is reduced by sea quark effects and the full 
QCD result becomes consistent with experiment. 

The results obtained with one-loop renormalization factors are subject to higher order corrections. However, 
as shown in Table IXXI VI for the quenched results, the difference between the perturbative and non-perturbative 
matchings is not large. This is because the 0{a) correction to the improved current in Eq. I|70|) is not large, and 
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the non-perturbative values for Za and 6a are close to those in tadpole improved perturbation theory. We may 
therefore expect that the uncertainty due to the perturbative matching is small also in the full QCD results. 

It is expected that various systematic uncertainties, including the scaling violation, would partially cancel in 
the ratio fx/ f-n- This expectation is supported by the good agreement of the quenched results between the 
perturbative and non-perturbative matching, which suggests that higher order corrections to the renormalization 
factors almost cancel in the ratio. Therefore the ratio is useful to discuss sea quark effects. Comparison of this 
quantity shows that the full QCD values are significantly closer to the experimental value ~ 1.22 by about two 
standard deviations than in quenched QCD. 



VIII. QUARK MASSES 



We calculate the up-down and strange quark masses through the axial vector Ward identity ( AWI) . The bare 
quark mass at simulation points is obtained by 

where and C]^ are the amplitudes of {A™^^ {t)P{Qy) and (P(i)P(0)'f) with the doubly smeared source and 
the local sink operators. 

We then carry out the chiral fit of the PS meson mass as a function of the AWI bare quark mass. The fitting 
function is obtained from Eq. with the replacement of the VWI masses with the AWI ones. We also drop 
all monomial terms in the sea quark mass, since the PS meson mass vanishes in the chiral limit "^q^^i = even 
for non-zero sea quark masses . The adopted form is 

_ E)PS,AWI AWI I ^PS,AWI / AWI\2 I ^PS,AWI AWI AWI ('70\ 

A*q,val "-^v lA^q,val7 "-^sv A^q.soa A'q.vab \'°) 

where 

Mqji = r^{K,,,)m^"^\Kso.,K,...K,..), (79) 

M^^i = r„(i^,ea)m^wi, (80) 

AWI _ 1 /'mAWI , AWI \ ^ 
"'■q.val — 2 V q^val,! ~r '"q,val,2j J \°'') 

Wq^aU = <^'(if«oa;i^val,„i^vau). (82) 

Our data and fit are shown in Fig. 221 

We adopt the fit H78|l because it is consistent with that in terms of the VWI quark mass (Eq. Ij^Sfl 'l. However, 
a function with fewer terms 

(T^ \^ (l^ -T^ V \\^ nPS.AWI ,,AWI , ,-<PS.AWI ,,AWI ,,AWI ^oq\ 

(.f^ol-f^scaj "IPSl-'^sca, ^val,l, ^val,2jj — -Dv Mq,val + '-'sv Mq.soa Mq.val \°-^) 

also gives an acceptable x^/dof. We use this as an alternative fit in our estimation of the systematic error due 
to the choice of the fitting function (see below). Parameters of these two fits arc summarized in Table IXXVI 

The bare AWI masses of the up-down and strange quarks are fixed in a way analogous to that described in 
Sec. IV CI by using Eq. (|75|l instead of Eq. . The matching to the MS scheme is made at the scale /i = 2/a 
using the one-loop renormalization constant 46, 47, 48,^,^J with the tadpole improvement. We use ap{q'^^), 
ap{ql^) as the expansion parameter in the one-loop expression of Za and bA, while we set q* — 2/a for other 
coefficients for which q* is not known. The MS quark mass is evolved to /i — 2 GeV using the four-loop beta 
function 

The VWI quark mass may differ from the AWI one because of explicit violation of chiral symmetry at finite 
lattice spacings. The difference between the AWI and VWI masses, therefore, gives insight into the size of scaling 
violation in our results. This leads us to repeat the calculation of quark masses using the VWI definition. The 
bare quark mass is calculated from K^, K^d and K, in Tables IVlIl and IXVll The MS mass is obtained by the 
one- loop matching at /i = 2/a and the four-loop running to /i = 2 GcV. The resulting AWI and VWI masses are 
summarized in Table IXXVII 
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For a calculation of quark masses in quenched QCD, the chiral extrapolation is carried out using a quadratic 
fit form obtained from Eq. H78() by dropping the third term which represents the sea quark mass dependence 

TTi {J^ r^r ^2 _ nPS.AWI „AWI , /^PS.AWI / ,AWIn2 (oa\ 

Obtained parameters are listed in Table IXXVIII We use either one-loop or non-perturbative value in Refs. [ssl 
Is^jH^iSSj for the renormalization factors. Numerical results are summarized in Table IXX Villi 

In Figure El our full QCD results are compared with estimates by the CP-PACS [13,113 and the QCDSF- 
UKQCD Collaborations [B^l- We observe good agreement for the AWI masses among the three groups even 
at the finite lattice spacing of a"-'^ ^ 2 GeV. These results are consistent also with the CP-PACS's result 
extrapolated to the continuum limit. This suggests that various uncertainties, such as scaling violation and 
higher order corrections to renormalization factors, are likely to be canceled in the ratio defining the AWI mass 
(Eq. G3). 

On the other hand, such a cancellation is not expected in the VWI mass. Indeed there is a sizable difference 
between our AWI and VWI results. We also observe that our and the CP-PACS results of VWI ra^d show a 
large discrepancy of about 18 % (6 standard deviations). These observations suggest that the scaling violation 
in our results of the VWI masses is not small. 

In both full and quenched QCD, therefore, we quote the AWI masses as the central value. We adopt ii'-input 
for TOg. The difference between vris with K- and 0- inputs is treated as the systematic error due to the choice of 
the meson mass input to fix m^. 

Additional systematic errors due to the choice of the chiral fit form and the uncertainty of the measured value 
of rg are included in our final results in full QCD. These errors are estimated in a similar way to that described 
in Sec. IV Cl bv using Eqs. (gUJ), (|^ and as alternative fit forms. 

Adding all errors in quadrature, we obtain 

m^(2GeV) = 3.223 (iro^^) (85) 

mf^(2 GeV) = 84.5 (tj^^") (86) 

ms/mud = 26.13 (i^;65) (87) 

in two-flavor QCD and 

m^(2 GeV) = 4.020 (0.077) (88) 

mp(2GeV) = 104.1(^24^1) (89) 

ms/mud = 2bm{tlf^). (90) 

in quenched QCD using the one-loop matching. The scaling violation is expected to be small and hence is 
ignored here. This point, however, should be checked in future studies. 

There is an additional uncertainty arising from the use of the perturbative value for ca, Za.p, and bA,p- 
Comparing the quenched result of Eqs. (|88|l - H90|) with those obtained by the non-perturbative matching given 

by 

m^(2 GeV) = 3.522(0.66) MeV (91) 

mf^(2 GeV) = 91.9 MeV (92) 

m,/m„<i = 26.08 , (93) 

we observe a systematic error of about 13 %. In quenched QCD the non-perturbative estimate of Z^/Zp = 1.19 
is very close to that in one-loop perturbation theory 1.22, since higher order corrections in Za and Zp partially 
cancel with each other. The non-perturbative value 6^ — 6p = 0.171 deviates significantly from that at one-loop 
—0.011. The 0{amq) term, however, is a small correction in our data. Most of the 13 % difference originates 
from the large deviation between the non-perturbative value ca = —0.083 in Ref. |33| and its one-loop value 
—0.013. Therefore, a non-perturbative determination of ca in full QCD is an important task toward a more 
precise calculation of the quark masses in future studies. 

In Fig. 1451 we compare the quark masses in full and quenched QCD. The chief observation is that sea quark 
effects reduce the light and strange quark masses by about 20 %. The magnitude of the sea quark effect is 
roughly consistent with the CP-PACS observation in Refs. p^lTH. 
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The quark mass ratio ms/mud in full QCD is consistent with the quenched value, because the sea quark 
effects in m„d and TOs almost cancel with each other in the ratio. We note that ms/mud — 26 is in a good 
agreement with the estimate of one-loop chiral perturbation theory 24.4(1.5) [Glj . 

Another important observation is that the deviation in between K- and (/>-inputs is reduced by effects 
of sea quarks: the deviation is 21 % (~ 24 MeV) in quenched QCD and 13 % (~ 12 MeV) in full QCD. This 
reflects the closer agreement of the meson spectrum in full QCD with experiment. The remaining deviation 
may be attributed to quenching of strange quarks and scaling violation. 



IX. CONCLUSIONS 



In this paper we have presented a high statistics study of the hadron spectrum and quark masses in two-flavor 
QCD using the plaquette gauge action and the fully 0(a)-improved Wilson quark action. We find firm evidence 
of sea quark effects at the simulated quark masses: the slopes dmy/dmpg and dmdec/dmpg are larger than 
in quenched QCD and the J parameter increases for lighter sea quarks. These findings do not suffer from 
systematic errors arising from the chiral extrapolation with respect to the sea quark mass, which is a major 
uncertainty particularly in recent studies with the Wilson-type quark action. Note that the use of a volume 
La> 1.8 fm at smaller sea quark masses TOps,sca/wv,soa — 0.6-0.7 is an important factor to control finite-size 
errors and in reaching our observations. 

The sea quark effect observed at the simulated quark masses means that the strange meson and baryon masses 
in full QCD show a better agreement with experiment than in quenched QCD. A similar reduction of quenching 
errors is also observed in the ratio Jk/ f-n- We also find that the sea quark effects lead to about 20 % reduction 
of quark masses. 

For baryons finite-size effects are large for the volume we used, which render sea quark effects unclear for 
lighter baryons. Further investigations on larger spatial volumes, of the order of 3 fm at the lightest sea quark 
mass, are needed to observe sea quark effects in the light baryons. 

The present work is carried out at a single lattice spacing of = 2.221(28) GeV. The 0(a)-improved Wilson 
quark action we employed is designed to have reduced scaling violation, and experiences in quenched QCD [62j 
support this expectation. Nonetheless scaling study of both the hadron spectrum and quark masses with this 
action is needed to establish the sea quark effects in the spectral quantities on a quantitative basis. 

Another important subject in future is simulations at much lighter sea quark mass, in particular below the 
p — > TTTT threshold. Such simulations will lead to a better control of the chiral extrapolation. This would also 
give insight into the chiral logarithmic singularity in the PS meson mass and decay constant, which we have 
not observed in our data. 
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APPENDIX A: HADRON MASSES 

Measured hadron masses are summarized in Tables IXXIXlHXXXI VI for full QCD, and Tables IXXXVMxn for 
quenched QCD. Our choice of the fitting range and resulting value of x^/dof are also shown in these tables. 
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TABLE I: Run parameters in simulations of two-fiavor QCD. The stop size At is given by the inverse of the number of 
the molecular dynamics steps (#MD). We denote the tolerance parameter in the stopping condition for the quark matrix 
inversion in calculations of the force and Hamiltonian by Af and Ah, respectively. CPU time required per trajectory 
on the full machine is written in units of minute. Number of measurement is denoted by Nmeas and the number of the 
exceptional configurations is written in the bracket. 



12^x48, A-HMC 


-^sca 


0.1340 


0.1343 


0.1346 


0.1350 


0.1355 


#MD 


68 


72 


80 


110 


175 


accept. 


0.667(8) 


0.667(10) 


0.668(10) 


0.772(6) 


0.799(9) 


Af 


10"* 


10"* 


10"* 


10"* 


10"* 


Ah 


10-" 


10-14 


10-14 


10-14 


10-14 




64.1(2) 


74.5(2) 


88.6(3) 


117.6(5) 


203.0(1.0) 


time/traj 


0.291 


0.351 


0.465 


0.829 


2.20 




3000 


3000 


3000 


3000 


3000 


^ meas 


300(2) 


300(0) 


300(1) 


300(2) 


300(7) 


mps/mv 


0.792(15) 


0.753(18) 


0.749(14) 


0.705(24) 


0.608(63) 


l()-''xlcS, IIMC 




0.1340 


0.1343 


0.1346 


0.1350 


0.1355 


#MD 


160 


160 


200 


200 


320 


accept. 


0.799(8) 


0.744(10) 


0.804(10) 


0.702(9) 


0.688(10) 


Af 


10-18 


10-18 


10-18 


10-18 


10-18 


Ah 


10-20 


10-20 


10-20 


10-20 


10-20 




105.9(2) 


124.0(2) 


148.0(5) 


203.0(7) 


362.2(1.9) 


time/traj 


1.92 


2.21 


3.86 


5.25 


14.7 


A^traj 


3000 


3000 


3000 


3000 


3000 


^ meas 


300(1) 


300(0) 


300(0) 


300(3) 


300(0) 


mps/mv 


0.802(5) 


0.779(12) 


0.752(8) 


0.707(13) 


0.586(21) 


2()-''xlcS, A-llMC 


-^^sca 


0.1340 


0.1343 


0.1346 


0.1350 


0.1355 


#MD 


100 


100 


115 


135 


250 


accept. 


0.673(7) 


0.627(8) 


0.670(8) 


0.663(9) 


0.755(8) 


Af 


10"* 


10"* 


10"* 


10"* 


10"* 


Ah 


10-14 


10-14 


10-14 


10-14 


10-14 




64.3(1) 


74.4(1) 


88.2(2) 


118.2(2) 


214.9(8) 


time/traj 


1.57 


1.79 


2.44 


3.75 


12.5 


^traj 


3000 


3000 


3000 


3000 


3000 


^ meas 


300(2) 


300(0) 


300(1) 


300(0) 


300(0) 


mps/mv 


0.802(5) 


0.773(5) 


0.745(18) 


0.705(5) 


0.586(8) 


2()-''xl8, S-li.\lC 


-^sca 


0.1340 


0.1343 


0.1346 


0.1350 


0.1355 


#MD 


80 


100 


100 


100 


160 


accept. 


0.676(5) 


0.771(5) 


0.749(4) 


0.666(6) 


0.678(7) 


Af 


10"® 


10"* 


10"* 


10"* 


10"* 


Ah 


10-14 


10-14 


10-14 


10-14 


10-14 




69.3(1) 


81.0(1) 


96.9(1) 


131.5(2) 


243.6(5) 


time/traj 


1.08 


1.53 


1.78 


2.35 


6.64 


^traj 


9000 


9000 


9000 


9000 


9000 


-^meas 


900(4) 


900(2) 


900(6) 


900(10) 


900(7) 


mps/mv 


0.799(1) 


0.779(1) 


0.753(2) 


0.709(2) 


0.600(4) 



22 



TABLE II: Parameters in Eq. 1191 and ro in full QCD. The first error is statistical. The second and third ones are the 
systematic error due to the choice of tmin and rmin. 



12^ X 48 




Vo 


a 


gc 


a 


ro 


0.1340 
0.1343 
0.1346 
0.1350 
0.1355 


0.549(19) (71) (37) 
0.552(14) (43) (47) 
0.568(14) (23) (34) 
0.603(10)(17)(13) 
0.658(10)(14)(29) 


0.301(13)(19)(21) 

0.288(9)(32)(44) 

0.288(10)(16)(64) 

0.296(8)(7)(19) 

0.333(7)(24)(21) 


0.76(5)(24)(12) 

0.74(5)(8)(11) 

0.67(4)(5)(8) 

0.57(3)(7)(4) 

0.47(3)(1)(7) 


0.0922(19)(30)(16) 
0.0838(10) (49) (36) 
0.0772(15) (18) (32) 
0.0675(11)(14)(14) 
0.0514(11)(25)(25) 


3.826(24) (35) (27) 
4.031(15)(76)(33) 
4.200(29) (23) (42) 
4.481(27)(59)(17) 
5.059(45) (80) (77) 


16^ X 48 


sea 


Vo 


a 




0" 


ro 


0.1340 
0.1343 
0.1346 
0.1350 
0.1355 


0.567(9) (18) (60) 
0.597(9) (5) (47) 
0.612(7) (23) (35) 
0.644(9) (8) (20) 
0.667(7) (10) (22) 


0.268(8)(3)(45) 

0.288(8)(13)(38) 

0.296(6)(21)(25) 

0.317(7)(5)(17) 

0.331(5)(17)(18) 


0.82(3)(9)(19) 
0.75(3)(2)(14) 
0.71(3)(5)(12) 
0.62(3)(7)(10) 
0.57(3)(1)(13) 


0.0919(10) (7) (44) 
0.0811(10)(18)(36) 
0.0737(8) (28) (22) 
0.0624(8)(9)(15) 
0.0515(6)(18)(10) 


3.880(14)(11)(31) 
4.098(16) (25) (34) 
4.287(19) (49) (24) 
4.621(21)(25)(26) 
5.059(23) (59) (32) 


20^ X 48 




Vo 


a 




a 


ro 


0.1340 
0.1343 
0.1346 
0.1350 
0.1355 


0.626(3)(12)(38) 
0.631(3)(3)(39) 
0.643(2)(21)(33) 
0.662(2) (4) (20) 
0.687(2)(6)(16) 


0.302(3)(10)(27) 

0.303(3)(1)(28) 

0.304(2)(15)(26) 

0.312(2)(6)(16) 

0.330(2)(8)(17) 


0.79(1)(3)(16) 

0.77(1)(1)(17) 

0.71(1)(7)(13) 

0.64(1)(1)(9) 

0.56(1)(1)(12) 


0.0866(4) (17) (25) 
0.0785(3) (6) (25) 
0.0716(3)(22)(21) 
0.0623(2)(16)(13) 
0.0509(2)(18)(13) 


3.946(5) (24) (17) 
4.143(6)(16)(23) 
4.336(6) (43) (24) 
4.635(7) (49) (20) 
5.092(8) (76) (33) 



TABLE III: Parameters in Eq. 11911 and ro in quenched QCD. 

lattice Vo a Qc u ro 

12=^ X 48 0.568(7)(17)(10) 0.276(5)(8)(13) 0.318(21)(36)(42) 0.0519(6)(17)(10) 5.149(24)(72)(22) 

16=* X 48 0.592(4)(12)(4) 0.278(3)(9)(11) 0.324(15)(29)(27) 0.0497(4)(11)(2) 5.255(15)(42)(15) 

20^ X 48 0.610(4)(5)(14) 0.287(3)(7)(18) 0.319(15)(9)(49) 0.0484(4) (2) (15) 5.309(16)(1)(47) 



TABLE IV: Autocorrelation time for plaquette (rp"™), PS meson propagator (rp™) and Wilson loop (tw™) for A-HMC 
and S-HMC simulations. All numbers are written in units of HMC trajectory. 



^^sca 


0.1340 


0.1343 


0.1346 


0.1350 


0.1355 


A-HMC 


cum 
' plaq 


11.4(1.5) 


35.4(9.8) 


21.9(4.7) 


10.3(2.4) 


13.6(2.2) 


cum 
T^PS 


13.7(2.6) 


13.5(3.1) 


17.6(4.3) 


7.7(0.9) 


8.5(3.1) 


cum 


5.0(0.4) 


5.9(0.8) 


5.9(1.1) 


4.6(0.7) 


6.2(0.9) 


S-HMC 


^^sea 


0.1340 


0.1343 


0.1346 


0.1350 


0.1355 


cum 
' plaq 


32.6(6.1) 


30.7(4.8) 


20.0(3.1) 


18.9(2.7) 


12.0(1.3) 


cum 


16.4(2.0) 


16.1(2.2) 


13.8(2.0) 


14.6(6.2) 


10.4(1.3) 


cum 
' W 


6.3(2.5) 


4.9(0.5) 


5.1(0.3) 


5.2(0.5) 


5.2(0.3) 



TABLE V: Parameters of diagonal fits to PS meson masses in method-A. We put "-" in some columns in this and the 
following tables when the corresponding term is not included in the fit. 



xVdof 




oPS 


^PS 
'-'diaK 




0.01 


0.136026(11) 


25.97(63) 


4.5(2.4) 




0.01 


0.136025(43) 


26.0(4.0) 


4(31) 


2(72) 
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TABLE VI: Parameters of diagonal fits to vector meson masses in method-A. 

X^/dof ^diaK -^diaR C'diag -Pdiag 

0.06 1.867(36) 0.248(17) -0.0083(19) 

0.04 1.90(10) 0.224(81) -0.002(21) -0.0005(17) 



TABLE VIL Parameters of combined chiral fits to PS meson masses in terms of VWI quark mass. 



xVdof 




Br 








^PS 


^PS 


0.33 


0.1360187(95) 


9.64(51) 


16.694(93) 


7.5(1.8) 


2.39(13) 


-6.68(41) 




0.33 


0.1360185(95) 


9.62(51) 


16.751(95) 


7.6(1.8) 


1.61(21) 


-6.78(41) 


0.63(11) 



TABLE VIIL Parameters of combined chiral fits to vector meson masses. 



xVdof 








cy 






0.55 


1.914(24) 


0.0349(51) 


0.1895(45) 






-0.00554(96) 


0.23 


1.877(33) 


0.051(14) 


0.1953(55) 


-0.0028(16) 


-0.00159(41) 


-0.00389(65) 



TABLE IX: Parameters of partially quenched chiral fits to vector meson masses. 

Ksc. xVdof ^PQ i^^Q 

0.1340 0.10(5) 0.5306(23) 0.6361(53) 

0.1343 0.08(5) 0.5067(23) 0.6648(58) 

0.1346 0.04(3) 0.4835(27) 0.6940(84) 

0.1350 0.14(10) 0.4394(29) 0.786(10) 

0.1355 0.27(20) 0.3885(34) 0.906(14) 



TABLE X: Parameters of combined chiral fits to octet baryon masses. 



xVdof 


AO 


B? 


F° 

V 


D° 


c? 


























0.82 


2.512(50) 


0.018(23) 


0.1286(32) 


-0.0344(34) 


0.0001(28) 


-0.00621(39) 


-0.00148(54) 






0.00047(23) 


0.00374(26) 


0.00057(35) 


0.00014(39) 










0.66 


2.489(51) 


0.020(23) 


0.1315(34) 


-0.0342(34) 


-0.0008(28) 


-0.00520(32) 


-0.00088(54) 


-0.00356(34) 




0.00074(25) 


0.00395(27) 


0.00044(35) 


0.00006(39) 











T.VIjLE XI: Paraiiictcrs of (■(niil)iii('(l cliiral fits to dccuplct l)aryc)ii iiiasst^s. 



xVdof 


A° 


BF 


B° 


CF 


CF 




0.92 


3.317(65) 


0.038(15) 


0.268(13) 






-0.0055(29) 


0.92 


3.327(85) 


0.049(36) 


0.253(16) 


-0.0012(43) 


0.0020(13) 


-0.0059(19) 



TABLE XII: Parameters of chiral extrapolation of Tq ^ . 



xVdof A 



1.5 -4.867(35) 0.6861(48) 



TABLE XIII: Parameters of diagonal fits to PS meson masses in method-B. 



x7dof 


K' 

c.diag 


D/PS 


t^diaK 




2.5 


0.136020(10) 


4.67(11) 


39.1(1.3) 




1.9 


().i:io<)7()(2:!) 


.-,.17(12) 


15(l;i) 


252(112) 



24 



TABLE XIV: Parameters of diagonal fits to vector meson masses in method-B. 

/Hof 4'^ R'^ r'^ n'^ 

1.4 0.3523(45) 1.481(38) -1.000(73) 

0.1 0.3379(89) 1.73(13) -2.19(59) 1.71(80) 



TABLE XV: Parameters of chiral extrapolation of ro in terms of PS meson mass squared. 

xVdof a;^ Bl^ Oo D'r„ 

0.01 5.629(34) -8.41(56) 16.3(2.7) -15.2(3.8) 



TABLE XVI: Hopping parameters corresponding to the light (Kud) and the strange quark mass with K- [Ks(K)) and 
0- input [Ks[(j>)) in full QCD. The lattice cutoff determined from p meson mass is also written. Error is statistical only. 

Kud Ks{K) Ks{<t>) [GeV] 

0.1359896(90) 0.134857(29) 0.134711(44) 2.221(28) 







TABLE XVII: Fit 


parameters of chiral extrapolations 


in quenched QCD. 




mis 


xVdof 
0.18 


Kc.q 

0.135315(12) 


2.859(35) 


3.09(31) 






niY 


xVdof 
0.04 


0.4146(54) 


0.797(22) 








rrioct 


xVdof 
0.66 


0.538(14) 


0.521(47) 


-0.158(44) 


-0.36(17) -0.059(92) 


Cq°'^ 
0.255(95) 


m^cc 


xVdof 
0.03 


^q 

0.712(14) 


B? 

1.123(59) 









TABLE XVIII: Hopping parameters corresponding to the light (Kud) and the strange quark mass with K- {Ks{K)) 
and input {Ks{(j})) in quenched QCD. 

Kud KsjK) Ks{4>) [GeV] 

0.135245(11) 0.133571(43) 0.133214(76) 1.834(24) 



TABLE XIX: J parameter calculated by Eq. I|68|l . We obtain J = 0. 4242(61) at A'sea = Jfud with the alternative definition 
l|69ll . 





0.1340 


0.1343 


0.1346 


0.1350 


0.1355 


Kud 


quenched 


J 


0.3827(22) 


0.3818(21) 


0.3801(34) 


0.3932(32) 


0.4018(37) 


0.4153(62) 


0.3735(66) 
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TABLE XX: Strange vector meson masses in GeV units in full and quenched QCD. The meson mass input to fix the 
strange quark mass is written in tlie brackets in the first column. The first error is statistical. The second and third 
ones are systematic error due to the choice of chiral fit forms and systematic uncertainty of ro- The deviation from the 
experimental spectrum is denoted by Am. 





m(iV/ = 2) 


Am{Nf = 2) 


m{Nf = 0) 


Am(A^/ = 0) 




0.8791(14)(t?)(i;^) 


-1.9% 


0.8706(16) 


-2.9% 


HK) 


0.9899(28)(l|«)(t^l) 


-2.9% 


0.9727(32) 


-4.7% 


K{4>) 


0.5272(31)(^f)(t^3?) 


+5.8% 


0.5466(39) 


+9.4% 




0.89390(0)(+173)(0) 


-0.2% 


0.89390(0) 


-0.2% 



TABLE XXI: Baryon spectrum in GeV units in full and quenched QCD. 





m{Nf = 2) 


Am(Af/=2) 


m(Nf = 0) 


Am(iV/ = 0) 


N 


1.015(19)(t^«)(l^») 


+7.7% 


1.004(26) 


+6.7% 


E{K) 


1.185{17){tll)itt) 


-0.7% 


1.162(20) 


-2.6% 


AiK) 


l.U3{18){tf){tl) 


+2.4% 


1.120(20) 


+0.4% 


E{K) 


1.291(16)(+14)(il) 


-1.8% 


1.260(14) 


-4.3% 


E(^) 


1.204(17)(+18)(t^) 


+1.0% 


1.194(19) 


+0.1% 


A(0) 


1.159(17)(+18)(1I) 


+3.8% 


1.146(19) 


+2.7% 


H(0) 


1.325(16)(+18)(t?) 


+0.7% 


1.315(11) 


+0.02% 


A 


1.328(25) (+24) (1^1) 


+7.5% 


1.318(26) 


+6.7% 


Y:*{K) 


1.433(21)(tf )(1^°) 


+3.5% 


1.414(22) 


+2.2% 


E*{K) 


1.537(17)(lf)(1^6) 


+0.4% 


1.510(18) 


-1.4% 


fl{K) 


1.642(14)(li«)(l«) 


-1.8% 


1.606(14) 


-4.1% 


E*(0) 


1.447(20)(ir)(l5) 


+4.4% 


1.436(20) 


+3.7% 


H*(0) 


1.566(15)(ir)(l|) 


+2.2% 


1.554(14) 


+1.4% 


n(0) 


1.684(ll)(lf )(1?) 


+0.7% 


1.671(10) 


-0.1% 



TABLE XXIL Parameters of combined chiral fits to PS meson decay constant in full QCD. 



xVdof 


A' 


Bl 


Bi 


cl 


a 


Csv 


0.001 


0.3415(71) 


0.0075(18) 


0.0349(18) 






-0.00128(42) 


0.001 


0.344(11) 


0.0019(55) 


0.0398(18) 


0.00046(70) 


-0.00105(19) 


-0.00069(23) 



TABLE XXIII: Parameters of chiral extrapolation to PS meson decay constant in quenched QCD. 



xVdof 






0.06 


0.07221(94) 


0.1246(42) 
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TABLE XXIV: Pseudoscalar decay constants in full and quenched QCD calculated with one-loop and non-perturbative 
(NP) matchings. Decay constants axe written in GeV units. The first error is statistical. The second and third ones for 
full QCD results are systematic error due to the choice of chiral fit forms and uncertainty of measured value of ro . 





one- loop ca, Za and b 


A 


NP Ca, Za and 6a 




Nf=2 


Nf=0 


Nf=0 




0.1372(31)(+25) (^") 


0.1337(24) 


0.1287(23) 


fK{K) 


0.1576(23)(+28) 


0.1497(20) 


0.1438(19) 




0.1603(21)(+29) (+3) 


0.1533(17) 


0.1473(17) 


fK{K)/U 


1.148(11) {tl') {tl) 


1.1195(63) 


1.1174(63) 


fK{ct>)/fn 


1.168(13) {tl') {It) 


1.1468(90) 


1.1443(89) 



TABLE XXV: Parameters of combined chiral fits to PS meson masses in terms of AWI quark mass in full QCD. 



xVdof 


n>PS,AWI 


^PS,AWI 




0.21 


18.08(12) 


0.33(19) 


-3.01(38) 


0.22 


18.13(11) 




-2.89(39) 



TABLE XXVI: Quark masses in two-flavor QCD in MS scheme at /i = 2 GeV. Values of m^d and rus are in MeV units. 
The meson mass input to fix rris is written in brackets. The quoted errors are statistical only. 





mud ms{K) ms{(j>) ms{K)/mud ms{4i)/mud 


AWI 
VWI 


3.223(43) 84.5(1.1) 96.4(2.2) 26.133(24) 29.78(38) 
2.488(57) 98.2(1.1) 110.4(2.3) 39.36(73) 44.21(93) 


TABLE XXVII: Parameters of chiral fit to PS meson masses in terms of AWI quark mass in quenched QCD. 


xVdof 


dPS.AWI ,^ps,awi 


0.22 




2.974(44) 1.31(61) 


TABLE XXVIII 
ruud and nis are 


Quark masses in quenched QCD in MS scheme at /i = 2 GeV calculated in quenched QCD. Values of 
in MeV units. Choice of input to fix ma is written in brackets. 




mud ms{K) msi^)) ms(K)/mud ms{(j))/mud 




one-loop matching 


AWI 
VWI 


4.020(77) 104.1(1.6) 128.2(3.3) 25.90(13) 31.88(56) 
4.628(77) 114.2(1.6) 136.9(3.2) 24.673(88) 29.57(42) 




non-perturbative matching 


AWI 
VWI 


3.522(66) 91.9(1.4) 113.1(2.9) 26.08(13) 32.11(56) 
4.315(72) 106.5(1.5) 127.7(3.0) 24.681(87) 29.58(42) 
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TABLE XXIX: Meson masses on 12^ x 48 lattice in iV/ = 2 full QCD. 



^ sea 


-ft^val.l 


-ft^val,2 


\tn 


in,tmaxl X"^/dof 




\tm 


in.imaxl X'^/dof 


mv 


U.lo4U 





1340 





1340 




9,24 


2.7 


0.619(10) 






9,24 


2.1 


0.782(12) 


0.1340 





1343 





1343 




9,24 


2.9 


0.592( 








9,24' 


2.3 


0.762 


13) 


0.io4U 





1346 





1346 




9,24 


2.9 


0.564( 


12 






9,24' 


2.5 


0.742 


13) 


0.1340 





1350 





1350 




9,24 


2.6 


0.525(12) 






9,24' 


2.7 


0.714 


1^ 


0.io4U 





1355 





1355 




9,24 


2.1 


0.4709(95 


) 




9,24' 


2.6 


0.679 


14) 


0.1340 





1358 





1358 




9,24 


1.9 


0.4344(79) 




9,24' 


2.4 


0.660 


13) 


0.1340 





1340 





1343 




9,24 


2.8 


0.606f 








9,24' 


2.2 


0.772 


13) 


0.io4U 





1340 





1346 




9,24 


3.0 


0.592( 


11) 






9,24' 


2.3 


0.762 


13) 


0.1340 





1340 





1350 




9,24 


3.0 


0.574^ 


121 






9,24' 


2.4 


0.749 


13) 


n 1 Q /I n 
U.io4U 





1340 





1355 




9,24 


2.7 


0.550( 


12) 






9,24 


2.5 


0.731 


13) 


0.io4U 





1340 





1358 




9,24 


2.4 


0.535( 


11) 






9,24 


2.5 


0.720 


14) 


U.lo4o 





1340 





1340 




9,24 


2.2 


0.5766 


/ ^ J 
^1 


) 




9,24 


2.6 


0.748 


13) 


O.io43 





1343 





1343 




'9,24 


2.2 


0.5474 


(51 


) 




9,24' 


2.7 


0.728 


14) 


O.io43 





1346 





1346 




9,24 


2.1 


0.5174 




) 




9,24' 


2.6 


0.704 


16) 


O.io43 





1350 





1350 




9,24 


2.2 


0.4755 


(50 






9,24' 


1.8 


0.665 


19) 


0.1343 





1355 





1355 




9,24 


2.4 


0.4173 


61 






9,24' 


0.9 


0.626 


16) 


O.io43 





1358 





1358 




'9,24 


2.9 


0.3798 


'69 






9,24' 


1.3 


0.613 


19) 


U.io4o 





1343 





1340 




'9,24 


2.2 


0.5621 


(53 






9,24' 


2.7 


0.738 


13) 


O.io43 





1343 





1346 




9,24 


2.1 


0.5325 


(49 






9,24' 


2.7 


0.716 




O.io43 





1343 





1350 




'9,24 


2.2 


0.5124 


(47 


) 




9,24' 


2.5 


0.699 


lei 


0.1343 





1343 





1355 




9,24 


2.2 


0.4867 




) 




9,24 


1.9 


0.675 


18) 


U.io4o 





1343 





1358 




9,24 


2.2 


0.4704 


(56 


) 




9,24 


1.3 


0.661 


18) 


U.lo4D 





1340 





1340 




9,24 


4.0 


0.5595 




) 




9,24 


2.5 


0.7122 




0.io4D 





1343 





1343 




'9,24 


4.7 


0.5304 


(49 






9,24' 


2.6 


0.6917 




0.1346 





1346 





1346 




9,24 


4.5 


0.5011 


(70 






9,24' 


2.6 


0.669 


11) 


0.1340 





1350 





1350 




'9,24 


3.0 


0.4596 


^80 






9,24' 


2.8 


0.636 


15) 
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TABLE XXXI: Meson masses on 16^ x 48 lattice in iV/ = 2 full QCD. 
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TABLE XXXII: Baryon masses on 16^ x 48 lattice in AT/ = 2 full QCD. 





ATval,! 


Aval. 2 


Aval, 3 


ftmin,imaxl X^/dof 




mE 


ftmin,tmaxl X^/dof 


rriA 


[imin,tmaxl X^ /dof 




nid 


2C 


0.1340 





1340 





1340 





1340 




10,24 


2.2 


1 


166 
















10,24 


5.1 


1 


256 




0.1340 





1343 





1343 





1343 




10,24 


2.3 


1 


128 


i!; 














10,24 


3.8 


1 


224 


??i 


0.1340 





1346 





1346 





1346 




10,24 


2.2 


1 


089 


20; 














10,24 


2.8 


1 


192 


15) 


0.1340 





1350 





1350 





1350 




10,24 


1.7 


1 


037 


18 














10,24 


2.0 


1 


149 


13) 


0.1340 





1355 





1355 





1355 




10,24 


1.3 





979 


18 














10,24 


1.4 


1 


096 


13) 


0.1340 





1358 





1358 





1358 




10,24 


1.7 





955 


18; 














10,24 


1.2 


1 


067 


16) 


0.1340 





1340 





1340 





1343 




10,24 


2.1 


1 


152 


20; 






10,24 


2.4 


1.155 


22; 






10,24 


4.9 


1 


242 


20) 


0.1340 





1340 





1340 





1346 




10,24 


2.0 


1 


138 


20 






10,24 


2.7 


1.144 


^23 






10,24 


4.3 


1 


231 


19) 


0.1340 





1340 





1340 





1350 




10,24 


1.9 


1 


120 


18, 






10,24 


3.0 


1.129 


24; 






10,24 


3.6 


1 


218 


17) 


0.1340 





1340 





1340 





1355 




10,24 


1.7 


1 


097 








10,24 


2.9 


1.109 


22' 






10,24 


2.5 


1 


202 




0.1340 





1340 





1340 





1358 




10,24 


1.5 


1 


083 


It: 






10,24 


2.3 


1.097 


18^ 






10,24 


1.9 


1 


193 


n 


0.1340 





1343 





1343 





1340 




10,24 


2.5 


1 


142 


22' 






10,24 


2.1 


1.139 








10,24 


4.4 


1 


231 




0.1340 





1346 





1346 





1340 




10,24 


2.8 


1 


118 


24^ 






10,24 


2.0 


1.113 








10,24 


3.5 


1 


211 




0.1340 





1350 





1350 





1340 




10,24 


2.9 


1 


083(24' 






10,24 


1.6 


1.077(17' 






10,24 


2.5 


1 


183(14) 



30 



TABLE XXXII: (Continvcd) 





-ft^val.l 




K Veil. 3 


1 t-min 5 








niA 


l^min , 






md 


2C 


0.1340 


0.1355 


0.1355 


0.1340 




10 


24 


2.2 


1.040(21 






10,24 


1.1 


1.035(16) 




10 


24 


1.5 


1.149(12) 


0.1340 


0.1358 


0.1358 


0.1340 




10 


24 


1.6 


1.019 


18' 






10,24 


1.2 


1.011(16) 




10 


24 


1.2 


1.130 




0.1343 


u 




u 


104U 


u 


104U 






24 


1 n 

-L .U 


1.100 
















1 

-LU 


24 


9 


n 
u 


1 
1 


ZUO 


lO) 


0.1343 


u 




u 


1040 


u 


1040 




1 n 




n 7 

U. ( 


1.063 


13' 














1 n 

lU 


94 
z^ 



z 


A 

4 


1 
1 


1 / D 


16' 


0.1343 





1346 





1346 





1346 




1 n 


24 


7 


1.025 


13^ 














1 

-LU 


24 


2 


9 


I 


150 


18* 


0.1343 


u 


lOOU 


n 
u 


lOOU 


u 


lOOU 




1 n 


24 


u.o 


0.972 


15^ 














1 n 

-LU 


24 


Q 



1 
1 


1 
1 


1 lU 


27' 


0.1343 


u 


lOOO 


n 

u 


lOOO 


u 


lOOO 




10 


24 


1.5 


0.902 


23' 














10 


24 


9 


p) 



1 
1 


044 

U44 


36* 


0.1343 


u 


lOOo 


n 

u 


lOOO 


n 

u 


lOOO 




10 


24 


2.0 


0.863 


30^ 














10 


24 


9 
z 


1 
1 


1 
1 


004 

UU4 


40) 

iu; 


0.1343 


u 


1 '^A'^ 
lo^o 


n 

u 


104:0 


n 

u 


1 '^4n 

104U 




10 


24 


0.8 


1.076 


13' 






10,24 


0.8 










10 


24 


9 
z 


Q 



1 
1 


1 R'^ 
100 


15) 


0.1343 


n 

u 


104:0 


n 

u 


1040 


u 


104U 




10 


24 


0.7 


1.050 


13^ 






10,24 


0.7 


1 n^ii 


11 






10 


24 


9 
z 


U 


1 
1 


lUU 


16) 


0.1343 


u 


1041:0 


n 
u 


1040 


n 
u 


1 "^pin 

lOOU 




10 


24 


0.6 


1.032 


13' 






10,24 


0.8 


1 .UoO 


12; 






10 


24 


9 
z 


Q 


1 
1 


1 Pi Pi 

100 


16' 


0.1343 


u 


1040 


n 
u 


1040 


n 

u 


1 "^PiPi 
lOOO 




10 


24 


0.8 


1.008 


14^ 






10,24 


1.0 


1.015 


13 






10 


24 


Q 



4 


1 
1 


1 '^Q 
loy 


19* 


0.1343 


u 


1040 


n 
u 


1040 


n 

u 


lOOo 




10 


24 


1.0 


0.993 


15^ 






10,24 


1.4 


1.003 


15. 






10 


24 


Q 



Pi 



1 
1 


1 9PI 

IZO 


23* 


0.1343 


n 
u 


1 '^/in 

lO^tU 


n 
u 


104U 


n 
u 


1040 




10 


24 


0.8 


1.087 


14' 






10,24 


0.9 


1.088 


13; 






10 


24 


9 
z 


1 
1 


1 
1 


1 Q9 
lyz 


15' 


0.1343 





104:U 


yj 


104U 




yJ 


1040 




10 


24 


0.7 


1.039 


13^ 






10,24 


0.7 


1.037 


13 






10 


24 


9 
z 


Q 



1 
1 


1 PIR 
100 


16* 


0.1343 


u 


lOOU 


n 

u 


lOOU 


n 
u 


1040 




10 


24 


0.9 


1.005 


13^ 






10,24 


0.7 


1.001 


14; 






10 


24 


Q 



4 


1 
1 


1 '^Pi 
100 


20* 


0.1343 


u 


1 Pi Pi 

1000 


n 
u 


1000 


n 
u 


1040 




10 


24 


1.4 


0.961 


17' 






10,24 


1.0 


0.953 








10 


24 





1 
1 


1 
1 


uyo 


32' 


0.1343 


n 


lOOo 


n 

u 


lOOO 


n 

u 


1040 




10 


24 


2.2 


0.932 


19^ 






10,24 


1.4 


0.925 








10 


24 


9 
z 


u 


1 
1 


uuo 


38* 


0.1346 


u 


104U 


u 


104U 


u 


104U 




1 1 

1 1 




9 4 


1.064 


10 














1 1 

L L 


94 
z^ 


1 
1 


Q 



1 
1 


1 piQ 
100 


16' 


0.1346 


u 


1040 


u 


1040 


u 


1 'iA 
1040 




11 


24 


9 


1.028 


12^ 














11 


24 


1 
1 


Q 



1 
1 


1 99 

izz 


16* 


0.1346 





101U 


yJ 


104U 





104U 




1 1 

1 1 


94 


9 9 


n QQi 

u. y y -L 


14' 














1 1 

L L 


94 
z^ 


]^ 


S 





0Q1 


17* 


0.1346 


u 


lOOU 


n 
u 


lOOU 


u 


lOOU 




1 1 

1 1 


94 


1 8 


0.938 


16^ 














1 1 

L L 


94 
z^ 


1 
1 


7 


1 
1 


UOO 


19 


0.1346 


n 
u 


1000 


n 
u 


1 Pi Pi 

1000 


n 
u 


1 "^Pipi 

1000 




11 


24 


1.1 


0.862 


20' 














11 


24 


1 
1 


9 
z 


1 
1 


00s 
UUo 


30' 


0.1346 


u 


lOOo 


n 
u 


lOOo 


n 
u 


lOOo 




11 


24 


1.7 


0.821 


30^ 














11 


24 


n 
u 


Q 


n 
u 


yoi 


40* 


0.1346 


n 

u 


104U 


n 

u 


104U 


n 
u 


1 '^4n 

104U 




11 


24 


2.3 


1.017 


11^ 






11,24 


2.3 


1.014(13' 






11 


24 


1 
1 


Q 



1 
1 


110 
1 lU 


18* 


0.1346 


n 
U 


104D 


n 
u 


104D 


n 
u 


1040 




11 


24 


2.2 


1.004 


12' 






11,24 


2.2 


1.002 


13; 






11 


24 


1 
1 


Q 



1 
1 


OQQ 

uyy 


18* 


0.1346 



u 


104:U 


yJ 


104U 


n 


lOOU 




11 


24 


2.0 


0.972 


15^ 






11,24 


2.1 


0.975 


13 






11 


24 


1 
1 


Q 


1 
1 


07P1 
u / 


18* 


0.1346 


n 

u 


104U 


n 

u 


104U 


n 

u 


1 "^PiPi 
lOOO 




11 


24 


1.8 


0.945 


15^ 






11,24 


1.9 


0.953 


U] 






11 


24 


9 
z 


n 

u 


1 
1 


0(^1 

UUl 


20* 


0.1346 


u 


104D 


n 
u 


104D 


n 
u 


lOOo 




11 


24 


1.7 


0.926 


15' 






11,24 


1.6 


0.939 


16] 






11 


24 


9 
z 


1 
1 


1 
1 


OPiPi 
UOO 


22* 


0.1346 


u 


104U 


n 

u 


104U 


n 

u 


104U 




11 


24 


2.4 


1.039 


12^ 






11,24 


2.4 


1.041 


10 






11 


24 


1 

1 


Q 




1 
1 


1 '^0 

lOU 


17* 


0.1346 


n 

u 


1040 


n 

u 


1040 


n 

u 


104U 




11 


24 


2.3 


1.015 


13^ 






11,24 


2.3 


1.016 


12; 






11 


24 


1 

1 


Q 



1 
1 


1 OQ 
luy 


17* 


0.1346 


u 


lOOU 


n 
u 


1 Pin 

lOOU 


n 
u 


104D 




11 


24 


2.0 


0.957 


14' 






11,24 


1.9 


0.954 


16; 






11 


24 


1 

1 


Q 

y 


1 
1 


0fi9 

UDZ 


19* 


0.1346 


u 


lOOO 


n 
u 


lOOO 


n 

u 


104U 




11 


24 


1.7 


0.908 


18'' 






11,24 


1.4 


0.901 


16 






11 


24 


1 

1 


7 


1 
1 


0'^4 

U04 


23* 


0.1346 


n 

u 


lOOO 


n 
u 


lOOO 


n 

u 


104U 




11 


24 


1.1 


0.873 


21^ 






11,24 


1.4 


0.869 


^20^ 






11 


24 


1 

1 


p; 



1 
1 


01 Q 
uiy 


30* 


0.1350 


u 


1 '^/in 

104U 


u 


104U 


u 


104U 




1 1 

1 1 


94 


9 1 


u. yUfj 


11 














1 1 

L L 


94 
z^ 


1 
1 


4 


1 
1 


0*^1 
UOl 


14 


0.1350 


u 


1040 


u 


1040 


u 


1040 






94 


1 7 


0.925 


10^ 














1 1 

L L 


94 
z^ 


1 
1 





u 


QQPi 

yyo 


15* 


0.1350 





1346 





1346 


Q 


1346 




11 


24 


-L .U 


0.884 


ll'' 














11 


24 


I 


2 





958 


16* 


0.1350 


u 


lOOU 


n 
u 


1 Pin 

lOOU 


n 
u 


1 "^pin 

lOOU 






94 


1 7 


0.828 


12' 














1 1 

L L 


94 
z^ 


1 

1 





u 


yuo 


18* 


0.1350 





lOOO 




yJ 


lOOO 




yJ 


lOOO 




11 


24 


1.1 


0.759 


16^ 














11 


24 


1 

1 


1 

1 


n 
u 


ozu 


33* 


0.1350 


n 
u 


lOOO 


n 
u 


lOOo 


n 

u 


lOOo 




11 


24 


1.0 


0.715 


25^ 














11 


24 


n 

u 


Q 

y 


n 
u 


( 00 


49* 


0.1350 


n 
u 


lOOU 


n 
u 


1 Pin 

lOOU 


n 
u 


104U 




11 


24 


1.5 


0.879 


13' 






11,24 


1.4 


0.8708(92 






11 


24 


1 

1 


9 
z 


n 
u 


y4 / 


16* 


0.1350 


n 
u 


lOOU 


n 

u 


lOOU 


n 

u 


1040 




11 


24 


1.6 


0.864 


13^ 






11,24 


1.5 


0.8577(98 


i 




11 


24 


1 

1 


9 
z 


n 
u 


yo4 


16* 


0.1350 


n 
u 


lOOU 


n 

u 


lOOU 


n 

u 


104U 




11 


24 


1.7 


0.849 


13^ 






11,24 


1.5 


0.845(11' 






11 


24 


1 

1 


Q 



n 

u 


Q90 
yzu 


17* 


0.1350 


n 
U 


lOOU 


n 
u 


1 Pin 

lOOU 


n 
u 


1 "^Pipi 

1000 




11 


24 


1.3 


0.802 


12' 






11,24 


1.7 


0.807 








11 


24 


1 

1 


9 
z 


n 
u 


/ D 


20* 


0.1350 


n 
u 


lOOU 


n 
u 


lOOU 


n 

u 


lOOO 




11 


24 


1.1 


0.787 


13^ 






11,24 


1.7 


0.792 


if 






11 


24 


1 

1 


n 

u 


n 
u 


00 i 


23* 


0.1350 


n 
u 


104U 


n 
u 


1 '^4n 

104U 


n 

u 


lOOU 




11 


24 


1.5 


0.916 


10^ 






11,24 


1.5 


0.9247(97) 




11 


24 


1 

1 


9 
z 


u 


yoo 


14* 


0.1350 


n 
u 


1040 


n 
u 


1040 


n 

u 


lOOU 




11 


24 


1.4 


0.890 


10' 






11,24 


1.5 


0.897 








11 


24 


1 

1 


9 
z 


n 
u 


yu4 


15* 


0.1350 


n 
u 


104U 


n 
u 


104U 


n 

u 


1 "^pin 

lOOU 




11 


24 


1.5 


0.864 


^11^ 






11,24 


1.6 


0.868 


I2', 






11 


24 


1 

1 


9 
z 


n 
u 


yoo 


16) 


0.1350 


u 


1000 


n 
u 


1 Pi Pi 

1000 


n 
u 


1 "^pin 

lOOU 




11 


24 


1.5 


0.784 


'17' 






11,24 


1.2 


0.778 








11 


24 


1 

1 


1 
1 


n 
u 


SPil 
oOi 


25' 


0.1350 


u 


lOOo 


n 
u 


lOOO 


n 
u 


lOOU 




11 


24 


1.3 


0.755 


17^ 






11,24 


1.2 


0.747 








11 


24 


n 

u 


c 




n 
u 


R^ 1 

oil 


36* 


0.1355 


u 


104U 


u 


104U 


u 


104U 




1 1 

1 1 


94 


9 7 

Z. f 


0.891 


16 














1 1 

L L 


94 
z^ 


1 
1 


1 
1 


u 


QS1 

yol 


17' 


0.1355 


u 


1040 


u 


1040 


u 


1040 




11 


24 


O. -L 


0.851 


18'' 














11 


24 


1 
1 


9 
Z 


u 


Q/1 S 

y4o 


22* 


0.1355 





1346 





1346 





1346 






94 


o.u 


n 81 4 


9n' 

> 














1 1 

L L 


94 
z^ 


1 


3 





919 


■^n* 

ou 


0.1355 





1350 





1350 





1350 




1 1 

L 1 


94 




0.771 


24 














1 1 

L L 


94 
z^ 


1 


4 





892 


39) 


0.1355 





1355 





1355 





1355 




11 


24 


2.0 


0.707 


^29; 














11 


24 


1 


9 





864 


53) 


0.1355 





1358 





1358 





1358 




11 


24 


1.4 


0.666 


47' 














11 


24 


2 


2 





806 


65) 


0.1355 





1355 





1355 





1340 




11 


24 


3.0 


0.791 


27 






11,24 


2.8 


0.761 


22' 






11 


24 


1 


4 





889 


39) 


0.1355 





1355 





1355 





1343 




11 


24 


2.8 


0.776 


26; 






11,24 


2.7 


0.749 


23^ 






11 


24 


1 


5 





883 


42) 


0.1355 





1355 





1355 





1346 




11 


24 


2.6 


0.760 








11,24 


2.6 


0.738 


24' 






11 


24 


1 


5 





877 


44) 


0.1355 





1355 





1355 





1350 




11 


24 


2.3 


0.737 


S: 






11,24 


2.4 


0.724 


27^ 






11 


24 


1 


6 





870 


47J 


0.1355 





1355 





1355 





1358 




11 


24 


1.8 


0.686 


30' 






11,24 


1.8 


0.703 


29; 






11 


24 


2 








845 


57) 


0.1355 





1340 





1340 





1355 




11 


24 


5.7 


0.818 


27 






11,24 


2.2 


0.839 


23 






11 


24 


1 


4 





923 


34) 


0.1355 





1343 





1343 





1355 




11 


24 


5.3 


0.794 


27; 






11,24 


2.6 


0.815 


25, 






11 


24 


1 


3 





907 


37) 


0.1355 





1346 





1346 





1355 




11 


24 


4.4 


0.773 


26; 






11,24 


2.8 


0.793 


27; 






11 


24 


1 


3 





893 


40) 


0.1355 





1350 





1350 





1355 




11 


24 


3.1 


0.746 


27 






11,24 


2.7 


0.759 


26 






11 


24 


1 


4 





878 


43) 


0.1355 





1358 





1358 





1355 




11 


24 


1.7 


0.692 


^32^ 






11,24 


1.6 


0.669 


^34^ 






11 


24 


2 


3 





825 


55) 



31 



TABLE XXXIII: Meson masses on 20^ x 48 lattice in iV/ = 2 full QCD. 
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TABLE XXXIV: Barvon masses on 20'' x 48 lattice in iVj- = 2 full QCD. 
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11 


24 


1.1 


0.8589 


34) 


0.1355 





1346 





1346 





1355 




11 


24 


1 


5 





7387 


27 






11,24 


1.2 


0.7485 


27 






11 


24 


1.1 


0.8341 


37) 


0.1355 





1350 





1350 





1355 




11 


24 


1 


4 





6990 


29 






11,24 


1.2 


0.7057 


29 






11 


24 


1.2 


0.8012 


44) 


0.1355 





1358 





1358 





1355 




11 


24 


1 


2 





6125 


48 






11,24 


1.2 


0.6073 


49* 






11 


24 


1.4 


0.7404 


89) 



33 



TABLE XXXV: Meson masses on 12^ x 48 lattice in quenched QCD. 



iVval,l 


-R'val.2 


[tmin,tmaxl X^/^of 


mps 


[tmin,tmaxl /dof 


mv 


U.iozoU 





13260 




10,24 


0.6 


0.4906(17 






10,24 


1.2 


0.6079(43) 


n 1 ooriA 

U.iozyu 





13290 




10,24 


0.6 


0.4616 








10,24 


1.3 


0.5875 




n 1 ooo 1 





13331 




10,24 


0.7 


0.4201 


20 






10,24 


1.3 


0.5606 




n 1 ooo /I 





13384 




10,24 


0.9 


0.3620 


24' 






10,24 


1.2 


0.5285(79) 


n 1 o /I oo 





13432 




10,24 


1.0 


0.3018 


30 






10,24 


0.9 


0.502(12) 


U.1>34d5 





13465 




10,24 


0.9 


0.2519 


38 






10,24 


0.4 


0.484(20 


rt 1 o o/^ A 





13290 




10,24 


0.6 


0.4763 


18 






10,24 


1.3 


0.5977 




rt 1 O A 





13331 




10,24 


0.7 


0.4562 


19 






10,24 


1.3 


0.5840 




U.IozdU 





13384 




10,24 


0.8 


0.4296 


20 






10,24 


1.2 


0.5669 


;56) 


A 1 OOiSA 





13432 




10,24 


0.9 


0.4048 


23 






10,24 


1.0 


0.5518 


65) 







13465 




10,24 


1.1 


0.3870 


26 






10,24 


0.8 


0.5396 


;79) 


0.13290 





13331 




10,24 


0.7 


0.4412 


19 






10,24 


1.3 


0.5740 


;52) 


0.13290 





13384 




10,24 


0.8 


0.4139 


2r 






10,24 


1.2 


0.5571 


60) 


0.13290 





13432 




10,24 


0.9 


0.3884 


24 






10,24 


1.1 


0.5422 


;70) 


0.13290 





13465 




10,24 


1.0 


0.3700 








10,24 


0.8 


0.5302 


f85) 


0.13331 





13384 




10,24 


0.8 


0.3919 


22^ 






10,24 


1.3 


0.5442 


67) 


0.13331 





13432 




10,24 


1.0 


0.3651 


25 






10,24 


1.1 


0.5297 


80) 


0.13331 





13465 




10,24 


1.0 


0.3458 


28 






10,24 


0.8 


0.5180 


'96) 


0.13384 





13432 




10,24 


1.0 


0.3332(27' 






10,24 


1.0 


0.5147(97) 


0.13384 





13465 




10,24 


1.0 


0.3121 


3o; 






10,24 


0.7 


0.504( 


12) 


0.13432 





13465 




10,24 


0.9 


0.2782 


33^ 






10,24 


0.6 


0.492( 


15 



TABLE XXXVL Baryon masses on 12^ x 48 lattice in quenched QCD. 



-ffval.l -?^val,2 


Ji'val,3 


l^min • 


tmaxl X^/dof 


rriT. 


ftmin,tmaxl X^'/dof 


mA 


l^min , 


tmaxl X^/dof 




"id 


DC 


0.13260 


13260 





13260 




11,24 


0.7 


0.9394 


(76 
















11 


24 


0.7 


1 


014 




0.13260 


13260 





13290 




11 


24 


0.7 


0.9264 


(79 


i 




11 


24 


0.7 


0.9294(78 


) 




11 


24 


0.7 


1 


003 


Hi 


0.13260 


13260 





13331 




11 


24 


0.8 


0.9083 


86 






11 


24 


0.8 


0.9159 


(83 






11 


24 


0.7 





989 


12) 


0.13260 


13260 





13384 




11 


24 


0.9 


0.8841 


98 


! 




11 


24 


0.9 


0.8983 


94 






11 


24 


0.7 





970 


13) 


0.13260 


13260 





13432 




11 


24 


1.1 


0.861(12) 




11 


24 


0.9 


0.882 


11) 




11 


24 


0.8 





952 


14) 


0.13260 


13260 





13465 




11 


24 


1.4 


0.841(14 




11 


24 


0.8 


0.869(13) 




11 


24 


1.0 





937 


16) 


0.13290 


13290 





13260 




11 


24 


0.8 


0.9180 


(81 






11 


24 


0.7 


0.9148(82) 




11 


24 


0.7 





993 


12) 


0.13290 


13290 





13290 




11 


24 


0.8 


0.9048 


85 


! 














11 


24 


0.7 





984 


12) 


0.13290 


13290 





13331 




11 


24 


0.9 


0.8865(92 


) 




11 


24 


0.9 


0.8911(90) 




11 


24 


0.7 





969 


13) 


0.13290 


13290 





13384 




11 


24 


1.0 


0.862( 








11 


24 


0.9 


0.873 


101 






11 


24 


0.7 





950 


14) 


0.13290 


13290 





13432 




11 


24 


1.2 


0.838( 








11 


24 


0.9 


0.857(12) 






11 


24 


0.7 





931 


15) 


0.13290 


13290 





13465 




11 


24 


1.5 


0.818(15) 






11 


24 


0.8 


0.843(14] 






11 


24 


0.9 





916 


18) 


0.133310 


13331 





13260 




11 


24 


0.9 


0.8884(92) 




11 


24 


0.9 


0.8804(95 


) 




11 


24 


0.7 





965 


13) 


0.133310 


13331 





13290 




11 


24 


0.9 


0.8750(96) 




11 


24 


0.9 


0.8703(98) 




11 


24 


0.6 





955 


13) 


0.133310 


13331 





13331 




11 


24 


1.0 


0.856( 


















11 


24 


0.6 





942 


14) 


0.133310 


13331 





13384 




11 


24 


1.1 


0.831( 


Si 






11 


24 


1.0 


0.838(12) 






11 


24 


0.6 





922 


15) 


0.133310 


13331 





13432 




11 


24 


1.2 


0.806( 


14 






11 


24 


1.0 


0.821 


14) 






11 


24 


0.7 





903 


17) 


0.133310 


13331 





13465 




11 


24 


1.6 


0.784( 


17 






11 


24 


0.8 


0.806 


16 






11 


24 


0.8 





887 


20) 


0.133840 


13384 





13260 




11 


24 


1.0 


0.849( 


12) 






11 


24 


1.1 


0.834 


13) 






11 


24 


0.6 





927 


15) 


0.133840 


13384 





13290 




11 


24 


1.0 


0.835( 


12 






11 


24 


1.1 


0.824 


13) 






11 


24 


0.6 





917 


16) 


0.133840 


13384 





13331 




11 


24 


1.0 


0.816( 


13 






11 


24 


1.1 


0.809 


13 






11 


24 


0.6 





903 


16) 


0.133840 


13384 





13384 




11 


24 


1.0 


0.789( 


14) 
















11 


24 


0.7 





887 


18) 


0.133840 


13384 





13432 




11 


24 


1.1 


0.762( 


17 






11 


24 


0.9 


0.771 


17) 






11 


24 


0.7 





866 


21) 


0.133840 


13384 





13465 




11 


24 


1.4 


0.736( 


21 






11 


24 


0.8 


0.755 


21 






11 


24 


0.8 





849 


25) 


0.13432 


13432 





13260 




11 


24 


0.9 


0.810( 


16) 






11 


24 


1.4 


0.788 


19) 






11 


24 


0.6 





890 


20) 


0.13432 


13432 





13290 




11 


24 


0.9 


0.795( 








11 


24 


1.3 


0.777 


19) 






11 


24 


0.6 





880 


21) 


0.13432 


13432 





13331 




11 


24 


0.9 


0.775( 


li! 






11 


24 


1.2 


0.761 


19 






11 


24 


0.6 





866 


22) 


0.13432 


13432 





13384 




11 


24 


0.9 


0.747( 


20 






11 


24 


1.0 


0.738(21 






11 


24 


0.7 





848 


23) 


0.13432 


13432 





13432 




11 


24 


0.8 


0.716( 


25 
















11 


24 


0.8 





832 


27) 


0.13432 


13432 





13465 




11 


24 


0.9 


0.682( 


33) 






11 


24 


0.7 


0.695(32) 






11 


24 


0.8 





810 


31) 


0.13465 


13465 





13260 




11 


24 


1.1 


0.774( 


22 






11 


24 


1.0 


0.746 


36) 






11 


24 


0.5 





858 


28) 


0.13465 


13465 





13290 




11 


24 


1.1 


0.758( 


23 






11 


24 


1.1 


0.733 


35 






11 


24 


0.6 





848 


28) 


0.13465 


13465 





13331 




11 


24 


1.0 


0.735( 


25) 






11 


24 


1.1 


0.713 


34) 






11 


24 


0.6 





833 


29) 


0.13465 


13465 





13384 




11 


24 


0.9 


0.703( 


30 






11 


24 


1.0 


0.683 


37) 






11 


24 


0.7 





813 


31) 


0.13465 


13465 





13432 




11 


24 


0.7 


0.665( 


40 






11 


24 


0.8 


0.653 


44 






11 


24 


0.8 





794 


34) 


0.13465 


13465 





13465 




11 


24 


0.6 


0.633( 


55) 
















11 


24 


1.0 





784 


38) 



34 



TABLE XXXVII: Meson masses on 16^ x 48 lattice in quenched QCD. 



-K-val.l 


-ft^val,2 


ftmin,tmaxl X'^/dof 


mps 


ftmin,tmaxl X^/dof 


niY 


U.iozbU 





13260 




10,24 


0.6 


0.48454(94) 




10,24 


1.3 


0.6024(20) 


n 1 oonn 

u.ioiyu 





13290 




10,24 


0.6 


0.45480 


(97) 




10,24 


1.3 


0.5805 




n 1 ooo 1 





13331 




10,24 


0.6 


0.4119(10' 






10,24 


1.2 


0.5506 


27) 


n 1 ooo ,1 
U.ioo84 





13384 




10,24 


0.6 


0.3511 


l^i 






10,24 


1.1 


0.5121 




n 1 o ,1 oo 





13432 




10,24 


0.9 


0.2871 








10,24 


0.9 


0.4778 


55) 


0.13465 





13465 




10,24 


1.3 


0.2332(23' 






10,24 


0.6 


0.4531 


88) 


O.IjzdU 





13290 




10,24 


0.6 


0.46982 


(95 






10,24 


1.3 


0.5915 


21) 


U.IozdU 





13331 




10,24 


0.5 


0.44919 


(99 


) 




10,24 


1.3 


0.5766 


23) 


U.lozoU 





13384 




10,24 


0.5 


0.4216 


10] 






10,24 


1.2 


0.5574 


27) 


U.ioZDU 





13432 




10,24 


0.6 


0.3956 


12 






10,24 


1.2 


0.5405 


32) 







13465 




10,24 


0.6 


0.3775 


14; 






10,24 


1.0 


0.5291 


38) 


0.13290 





13331 




10,24 


0.5 


0.4337 


10] 






10,24 


1.3 


0.5656 


25) 


0.13290 





13384 




10,24 


0.5 


0.4054 


11 






10,24 


1.2 


0.5464 


28) 


0.13290 





13432 




10,24 


0.6 


0.3786 


12; 






10,24 


1.1 


0.5294 


34) 


0.13290 





13465 




10,24 


0.7 


0.3597 


14] 






10,24 


1.0 


0.5180 


40) 


0.13331 





13384 




10,24 


0.6 


0.3825 


11 






10,24 


1.1 


0.5314 


31) 


0.13331 





13432 




10,24 


0.7 


0.3543 


12 






10,24 


1.1 


0.5143 


37) 


0.13331 





13465 




10,24 


0.7 


0.3343 


14^ 






10,24 


1.0 


0.5029 


45) 


0.13384 





13432 




10,24 


0.7 


0.3206(133 






10,24 


1.0 


0.4950(44) 


0.13384 





13465 




10,24 


0.8 


0.2986 


16; 






10,24 


0.9 


0.4835 




0.13432 





13465 




10,24 


0.9 


0.2625 


19^ 






10,24 


0.8 


0.4663 





TABLE XXXVIII: Baryon masses on 16^ x 48 lattice in quenched QCD. 



Kva].! i^val.2 


-^^val.a 


1 ^min ; 


tmaxl X'^/dof 




[tmin,imaxl X''/dof 


ruA 


[tmin,tmaxl X^/dof 


TTldcc 


0.132600 


13260 





13260 




11 


24 


0.9 


0.9081 


(44 














11.24 


1 


7 


0.9865 




0.132600 


13260 





13290 




11 


24 


0.9 


0.8948 


(45 






11.24 


0.9 


0.8969(45 


) 




11.24 


1 


7 


0.9746 




0.132600 


13260 





13331 




11 


24 


0.9 


0.8762 


46 






11.24 


0.9 


0.8817 


(47 






11.24 


1 


7 


0.9602 


59) 


0.132600 


13260 





13384 




11 


24 


0.8 


0.8512 


49 






11,24 


0.9 


0.8622 


51 






11.24 


1 


7 


0.9415 


64) 


0.132600 


13260 





13432 




11 


24 


0.7 


0.8273 


(56 






11,24 


0.9 


0.8454 


(56 






11,24 


1 


8 


0.9247 


(72) 


0.132600 


13260 





13465 




11 


24 


0.6 


0.8107 


67 






11,24 


0.8 


0.8351 


(65 






11.24 


1 


6 


0.9146 


86 


0.132900 


13290 





13260 




11 


24 


0.9 


0.8847 


46 






11,24 


0.9 


0.8824 


46 






11,24 


1 


7 


0.9642 


58 


0.132900 


13290 





13290 




11 


24 


0.9 


0.8712 


(47 














11,24 


1 


7 


0.9549 


(60) 


0.132900 


13290 





13331 




11 


24 


0.9 


0.8523 


48 






11,24 


0.9 


0.8559 


(49 






11,24 


1 


7 


0.9388 


63) 


0.132900 


13290 





13384 




11 


24 


0.8 


0.8266 


(52 






11,24 


0.9 


0.8363 


53 






11,24 


1 


7 


0.9197 


69) 


0.132900 


13290 





13432 




11 


24 


0.7 


0.8019 


(61 






11,24 


0.9 


0.8192 


(59 






11,24 


1 


6 


0.9025 


(78) 


0.132900 


13290 





13465 




11 


24 


0.6 


0.7849 


(73 






11,24 


0.8 


0.8086 


(69 






11,24 


1 


5 


0.8922 




0.133310 


13331 





13260 




11 


24 


0.9 


0.8526 


(50 






11,24 


0.9 


0.8461 


49 






11,24 


1 


7 


0.9351 




0.133310 


13331 





13290 




11 


24 


0.9 


0.8388 


(51 






11,24 


0.9 


0.8348(51 






11,24 


1 


7 


0.9242 


67) 


0.133310 


13331 





13331 




11 


24 


0.9 


0.8193 


(53 














11,24 


1 


6 


0.9108 


72) 


0.133310 


13331 





13384 




11 


24 


0.8 


0.7925 


(59 






11,24 


0.9 


0.7994(58 


) 




11,24 


1 


6 


0.8893 


78) 


0.133310 


13331 





13432 




11 


24 


0.7 


0.7665 


70 






11,24 


0.9 


0.7819 


(66 






11,24 


1 


4 


0.8715(91) 


0.133310 


13331 





13465 




11 


24 


0.7 


0.7490 


84 






11,24 


0.8 


0.7706 


78 






11,24 


1 


3 


0.861( 


}^\ 


0.133840 


13384 





13260 




11 


24 


0.9 


0.8111 


(58 






11,24 


0.7 


0.7954 


(60 






11,24 


1 


6 


0.8969 


78) 


0.133840 


13384 





13290 




11 


24 


0.9 


0.7967 


60 






11,24 


0.8 


0.7838 


(62 






11,24 


1 


5 


0.8855 


82) 


0.133840 


13384 





13331 




11 


24 


0.9 


0.7761 


64 






11,24 


0.8 


0.7680 


66 






11,24 


1 


4 


0.8697 


(88) 


0.133840 


13384 





13384 




11 


24 


0.9 


0.7475 


74 














11,24 


1 


3 


0.852( 


10) 


0.133840 


13384 





13432 




11 


24 


0.8 


0.7195(90 






11,24 


0.9 


0.7294(85) 




11,24 


1 


2 


0.831( 


12) 


0.133840 


13384 





13465 




11 


24 


0.8 


0.701( 


11) 




11,24 


0.9 


0.717(10) 




11,24 


1 


1 


0.821( 


15) 


0.134320 


13432 





13260 




11 


24 


1.0 


0.7749(75 


) 




11,24 


0.7 


0.7442(84 


) 




11,24 


1 


3 


0.862( 


11) 


0.134320 


13432 





13290 




11 


24 


1.0 


0.7596 


(79 






11,24 


0.7 


0.7324 


(88 






11,24 


1 


2 


0.851( 


11) 


0.134320 


13432 





13331 




11 


24 


1.0 


0.7376 


(86 






11,24 


0.7 


0.7161 


95 






11,24 


1 


1 


0.834( 


13) 


0.134320 


13432 





13384 




11 


24 


0.9 


0.707( 


10) 






11,24 


0.8 


0.695(11) 




11,24 


1 





0.814( 


15) 


0.134320 


13432 





13432 




11 


24 


0.9 


0.677( 


13) 














11,24 


1 





0.800( 


18) 


0.134320 


13432 





13465 




11 


24 


0.9 


0.656( 


16) 






11,24 


1.0 


0.665(16) 






11,24 


1 
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TABLE XXXIX: Meson masses on 20^ x 48 lattice in quenched QCD. 
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TABLE XL: Baryon masses on 20^ x 48 lattice in quenched QCD. 
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FIG. 1: Effective mass of PS meson at Jfsoa = ifvai = 0.1340 (left figure) and 0.1355 (right figure) on 20^ x 48 lattice in 
full QCD. We use the local sink operator for all data. 
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FIG. 2: Effective mass of vector meson at Ksea. = Kvai = 0-'i-S40 (left figure) and 0.1355 (right figure) on 20^ x 48 lattice 
in full QCD. We use the local sink operator for all data. 
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FIG. 3: Effective mass of octet baryon at iirsea = -?fvai = 0.1340 (left figure) and 0.1355 (right figure) on 20^ x 48 lattice 
in full QCD. We use the local sink operator for all data. 
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FIG. 4: Effective mass of decuplct baryon at Ksa, = K^^^i =0.1340 (left figure) and 0.1355 (right figure) on 20* x 48 
lattice in full QCD. We use the local sink operator for all data. 
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FIG. 5: Double exponential fits to PS meson (left figure) and octet baryon masses (right figure) at ii'sea = 0.1355. Right 
panel in each figure shows fitted masses determined from double exponential (open symbol) and single exponential fit 
(filled symbol). The local sink operator is used for all data. 



r = V5 



10 



1.02 



r = 5 



FIG. 6: Effective potential energies VeB{r,t) as a function of temporal separation t at ii'sea = 0.1350 on 20* x 48 lattice. 
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FIG. 7: Static quark potential on 20 x 48 lattice. Left and right figures show data at Jfsea = 0.1340 and 0.1355, 
respectively. 
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FIG. 8: Sommer scale on 20^ x 48 lattice as a function of rmin- Left and right figures show data at iC3ea = 0.1340 and 
0.1355 
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FIG. 9: Bin size dependence of jack-knife error of hadron masses on 20^ X 48 lattice in full QCD. 
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FIG. 10: Bin size dependence of jack-knife error of static potential at r = 5 on 20^ x 48 lattice in full QCD. 
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FIG. 11: Bin size dependence of jack-knife error of meson (left figures) and baryon masses (right figures) with Kvmi- 
0.13432, which corresponds to mps.vai/jTiv.vai — 0.6, on 20'^ x 48 lattice in quenched QCD. 
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FIG. 12: Bin size dependence of jack-knife error of static potential at r = 5 on 20 x 48 lattice in quenched QGD. 
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FIG. 13: Diagonal data of PS (left figure) and vector meson masses (right figure) as a function of 20^/V^. Meson 
masses in the infinite volume limit at i^sca = 0.1340 and 0.1350 are determined by the constant fit to data on two larger 
volumes, while we assume the linear dependence Eq. ()27^ at T^'soa = 0.1355. We also plot the prediction from the analytic 
formula |3ql for the PS meson mass at the lightest sea quark mass, by long dashed line. 
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FIG. 14: Diagonal data of octet (left figure) and decuplet baryon masses (right figure) as a function of 20^ /V. 
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FIG. 15: Volume size dependence of meson (left figure) and baryon masses (right figure) at i^vai = 0.13432 in quenched 
simulations. 




FIG. 16: Valence quark mass dependence of relative difference between hadron masses measured on 16^ and 20^ lattices. 
Left and right figure show data for the PS meson and the octet baryon, respectively. We define i^vai.ave by l/JS'vai,ave = 

(l/T^val.l + l/Ji'val,2)/2 for mCSOR masses m(A'sea; i^val.l, i^val,2) and l/JSTval.ave = (l/JS'val,! + l/.f!^val,2 + l/^fval,3)/3 for 

baryon masses m{Ksei,; ATvai,!, /4rvai,2, ^^vai.a). 



42 




FIG. 17: Vector meson (left figure) and decuplet baryon masses (right figure) at 1^863 = 0.1355 on three spatial volumes 
as a function of PS meson mass squared. Linear fit curve to each data set is shown as a guide for eyes. 




FIG. 18: Relative size of FSE for meson (left figure) and baryon masses (right figure) at /fsca = 0.1355 on 20^ x 48. Dot- 
dashed and dashed lines show the location where i^vai.ave = ^^'sea and 0.1350, respectively. The latter roughly corresponds 
to the strange quark mass. 
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FIG. 19: Magnitude of FSE on octet baryon mass at iC'soa^ 0.1355 as a function of spatial linear extent in lattice units. 
Circles are results for diagonal data, while squares represent those at ifvai = 0.1350. 




FIG. 20: Test of logarithmic singularity in quark mass dependence of PS meson mass. We use the quark mass defined 
through the axial vector Ward identity in this plot. Solid and dashed lines are fit curves of Eq. H29^ assuming / to be a 
free parameter or fixed to the experimental value. 




FIG. 21: Comparison of quadratic and cubic diagonal fits in method-A. 
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FIG. 22: Combined chiral extrapolation of PS meson masses in terms of VWI quark mass. 
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FIG. 23: Combined chiral extrapolation of vector meson masses. 
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FIG. 24: Combined chiral extrapolation of octet baryon masses. Top figure shows extrapolation of diagonal data. 
Middle and bottom figures are data at Ksea. = Kva,i,i = 0.1340 and 0.1355 for E-like (left panels) and A-like baryons 
(right panels). 
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FIG. 25: Combined chiral extrapolation of decuplet baryon masses. 
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FIG. 27: Comparison of quadratic and cubic diagonal fits in method-B. 
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FIG. 28: Relative magnitude of contribution of linear and quadratic terms in quadratic diagonal fit of PS meson masses. 
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FIG. 29: Relative magnitude of contribution of constant, linear and quadratic terms in quadratic diagonal fit of vector 
meson masses. 
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FIG. 30: Comparison of diagonal fits in method-A and B. 
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FIG. 31: Chiral extrapolation of ro in terms of mpg , 
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FIG. 32: Comparison of sea quark mass dependence of romv (left figure) and ro in lattice units (right figure) with those 
in CP-PACS data. We estimate romv in the chiral limit in the CP-PACS data from linear fit in terms of (romps)^. 




FIG. 33: Chiral extrapolation of meson (top figures) and baryon masses (bottom figures) in quenched QCD. For octet 
baryon masses, we plot only data at ifvai,2 = 0.13260 and 0.13432 for simplicity. 
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FIG. 34: Vector meson mass as a function of PS meson mass squared at each sea quark mass in full QCD and in quenched 
QCD. The experimental values of meson masses are also plotted using our result ro = 0.497 fm, which is determined from 
Eq. and Kud and a in Table IXVTI 
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FIG. 35: J parameter defined by Eq. I|68|l in full (left panel) and quenched QCD (right panel). Dashed lines are 
reproduced from combined chiral fit, Eq. 14111 . We also plot values calculated from an phenomenological definition, 
Eq. H69|l . using experimental spectrum (open diamond) and our results in Table IXXl (^filled square). 
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FIG. 36: Decuplet baryon masses as a function of PS meson mass squared at each sea quark mass in full and quenched 
QCD. 
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FIG. 37: Comparison of strange meson masses between full and quenched QCD. Experimental values are shown by 
horizontal lines. 
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FIG. 38: Mass of K* meson with if-input as a function of lattice spacing. 
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FIG. 39: Baryon spectrum with Jf-input. Experimental values are shown by horizontal lines. 
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FIG. 40: Nucleon mass (bottom panel) and H baryon mass with ii'-input (top panel) as a function of lattice spacing. 
Open triangles represent the CP-PACS results in quenched QCD using the standard plaquette gauge and the Wilson 
quark actions, while open squares are obtained with the renormalization group improved gauge and the tadpole improved 
clover actions. 




FIG. 41: Combined chiral extrapolation of PS meson decay constants. 
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Combined chiral extrapolation of PS meson masses in terms of AWI quark mass. 
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FIG. 44: Comparison of light (left figure) and strange quark mass with if-input (right figure) in two-flavor QCD. 
Triangles represent the CP-PACS results of the VWI quark mass using Kc determined by partially quenched chiral 
extrapolations. The CP-PACS result in the continuum limit was obtained by combined linear extrapolation of three 
data. We note that SESAM-T^L's results in Ref. are consistent with these results within large error arising from 
their continuum extrapolation. 
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FIG. 45; Comparison of light (left panel) and strange quark mass (right panel) in full and quenched QCD. 



